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ABSTRACT 



The feasibility of the breed/burn fuel cycle concept, 
originally proposed for Pu239/U FBR applications, has been 
evaluated for use in a U233/Th fueled pebble bed HTGR. Be- 
cause of its excellent neutron economy, its versatility with 
respect to fueling strategies and the high fuel burnup 
achievable as well as the ability to approach continuous 
fueling operation, the 'pebble bed HTGR w as judged to be the 
best suited, of reactors currently of commercial interest, 
for implementation of the breed/burn fuel cycle in a thermal/ 
epithermal system. 

Three necessary (but not sufficient) criteria for success- 
ful system operation were defined: (1) achieving system criti- 
cality, k 1.0; (2) achieving a breed zone fissile production 

rate such that the asymptotic enrichment attained in the breed 
zone ai.d the end-of-life burnup in the burn zone (expressed 
in percent) are related by e^gy > B EOL ; an ^ (3) accounting for 
zone-to-zone leakage by providing adequate excess reactivity, 
shown to require that k a 1.5. 

State-of-the-art computation methods, PDQ -7/AMPX, were 
employed for the analysis of homogeneous breed/burn pebble bed 
HTGR regions. Breed zones achieved their highest asymptotic 
enrichment (^3.5%) for carbon-to-heavy metal atom ratios, 

C/HM - 110; when charged to the burn zone this maximized 
k u , ( , v-1.3) . Highly-moderated burn zones (C/HM - 450) were most 
able to attain high fuel burnups (^100,000 MWD/MTHM) while 
sustaining the highest possible values of k a . . No combination 
of the compositions evaluated here (110 _< C/HM _> 450), however, 
were able to achieve the k ro (k^, ^>1.5) required to overcome the 
effect.: of zone- to- zone leakage. For these compositions then, 

successful operation of a breed/burn system was not possible. 
Further study of reactor systems operating with very hard breed 
zone eiptherma 1 spectra (C/HM < 110) should be undertaken if a 
feasible breed/burn system is to be pursued since it was found 
that asymptotic zone enrichment reaches a maximum in a spectrum 
between that typical of a hard-spectrum HTGR and soft-spectrum 
LMFBR. 
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ABSTRACT 



The feasibility of the breed/burn fuel cycle concept, 
originally proposed for Pu239/U FBR applications, has been 
evaluated for use in a U233/Th fueled pebble bed HTGR. Be- 
cause of its excellent neutron economy, its versatility with 
respect to fueling strategies and the high fuel burnup 
achievable as well as the ability to approach continuous 
fueling operation, the pebble bed HTGR was judged to be the 
best suited, of reactors currently of commercial interest, 
for implementation of the breed/burn fuel cycle in a thermal/ 
epithermal system. 

Three necessary (but not sufficient) criteria for success- 
ful system operation were defined: (1) achieving system criti- 
cality, k .>1.0; (2) achieving a breed zone fissile production 

rate such that the asymptotic enrichment attained in the breed 
zone and the end-of-life burnup in the burn zone (expressed 
in percent) are related by ®EOL» an< 3 (3) accounting for 

zone-to-zone leakage by providing adequate excess reactivity, 
shown to require that k ro >_ 1.5. 

State-of-the-art computation methods, PDQ -7/AMPX, were 
employed for the analysis of homogeneous breed/burn pebble bed 
HTGR regions. Breed zones achieved their highest asymptotic 
enrichment (v3.5%) for carbon-to-heavy metal atom ratios, 

C/HM = 110; when charged to the burn zone this maximized 
koo (vl.3) . Highly-moderated burn zones (C/HM ^ 450) were most 
able to attain high fuel burnups (^100,000 MWD/MTHM) while 
sustaining the highest possible values of koo . No combination 
of the compositions evaluated here (110 £ C/HM _> 450), however, 
were able to achieve the k^ (k^ ^>1.5) required to overcome the 
effects of zone-to-zone leakage. For these compositions then. 
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successful operation of a breed/burn system was not possible. 
Further study of reactor systems operating with very hard breed 
zone eipthermal spectra (C/HM < 110) should be undertaken if a 
feasible breed/burn system is to be pursued since it was found 
that asymptotic zone enrichment reaches a maximum in a spectrum 
between that typical of a hard-spectrum HTGR and soft-spectrum 
LMFBR. 
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I . INTRODUCTION 



1 . 1 Foreword 

In the past the general public has understood the term 
"nuclear proliferation" in a restricted sense, referring almost 
exclusively to the widespread deployment of nuclear weapons 
throughout the world. The realization that a growing number of 
foreign countries are operating commercial nuclear powered 
electrical generating stations and the implied relationship be- 
tween nuclear electric plants and weapons has resulted in a 
broader interpretation of proliferation. In the United States, 
this interpretation has achieved the status of official govern- 
ment policy and has led to both national and international 
reassessments, e.g. the Nonproliferation Alternative System 
Assessment Program (NASAP) and the International Fuel Cycle 
Evaluation (INFCE) efforts. The question of the safeguarding 
of sensitive nuclear materials has become a major area of 
emphasis in furthering the development and application of the 
nuclear fuel cycle for commercial electrical production. 

Previously it was envisioned that light water reactors 
(LWR) would serve an interim role as precursors to fast breeder 
reactors (FBR) . This status is testified to by the omission in 
most of the older long range systems studies of "advanced LWRs . " 
Until such time as FBRs could become commercially feasible and 
sufficient fissile reserves could be bred to permit their 
widespread deployment, LWRs were to be the mainstay of commer- 
cial nuclear power. 
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Currently, however, the transition to a breeder economy 
is not envisioned for the United States, at least in the 
short term, so it is prudent to consider other options to 
meet the needs for a growing electrical generating capacity 
[Nl] . Since the objectionable feature of the FBR is that it 
utilizes fertile uranium-238 (U238) to breed large amounts of 
fissile plutonium, new options must permit utilizing nuclear 
fuels in such a manner that the production of and/or commer- 
cial traffic in materials suitable for weapons applications is 
minimized at all levels. 

The dilemma posed by this objective is highlighted by 
rising oil prices and limited and unevenly distributed national 
energy resources. Many nonnuclear weapons states, desiring to 
provide for their own growing energy needs, have opted to 
develop civilian nuclear power programs. In the United States, 
the lack of expansion of enrichment capacity or the commercial 
deployment of reprocessing and mixed oxide fabrication capabil- 
ity has increased the economic and commercial incentives for 
other countries to develop such capabilities independently [PI] 

In view of the above considerations, development and evalu 
ation of the commercial feasibility of schemes utilizing the 
nuclear fuel cycle in a manner so as to minimize proliferation 
risks, while providing suitable economic and commercial incen- 
tives for their widespread deployment, has become an important 
national and international objective. It is the purpose of the 
present thesis to evaluate one particular class of prolifera- 
tion-resis tent fueling schemes; here designated the "breed/burn 
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concept. As will be discussed further, in its fullest develop- 
ment this scheme would allow charging only fertile material to 
a reactor after initial startup and dispatch of fully spent 
fuel to disposal rather than reprocessing and reuse. 
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1 . 2 Background 

Throughout its history, commercial nuclear power in the 
United States has been subjected to several periods of chang- 
ing emphasis. To the public utilities, short term fuel cycle 
costs have always been a most important consideration. While 
long-term trends are reviewed and evaluated, the utility 
always seeks to minimize the current costs of electricity to 
the consumer, hence causing the short-term economic consider- 
ations to dominate. Indeed, the very practice of present 
value discounting emphasizes the short-term. 

National energy planning, on the other hand, is concerned 
with being able to insure the availability of adequate energy 
resources within reasonable economic limits for the very long- 
term. With specific application to the nuclear fuel cycle, 
although most quantities of high-cost fissile uranium in the 
earth's crust can conceivably be considered a source of energy, 
the cost of utilizing it is prohibitive. In this light, a 
more economically viable solution to the long-term energy 
supply problem is to produce most future fissile materials by 
irradiating fertile materials in the fuel elements of a nuclear 
reactor. However, to accomplish this without relying on the 
commercialization of the FBR, it becomes necessary to operate 
other types of reactors as near- or net-breeders. 

Current pressurized water reactors (PWR) operate with 
conversion ratios on the order of 0.60 [Dl] on the uranium-238 
fuel cycle. High temperature gas-cooled reactors (HTGR) , how- 
ever, typically operate with thorium-232 (Th232) to uranium-233 



-15- 



(U233) conversion ratios of 0.76 [Dl] , based on optimization 
of near term economics. Conversion ratios of 0.84 are feasible 
and favored for ore conservation, however [Dl] , and even 
higher values are readily obtainable [L4,R2,S3]. Thus, we will 
focus on this inherently superior class of reactors in the pre- 
sent study. Most HTGRs built or designed to date utilize a 
highly enriched uranium/ thorium (HEU/Th) fuel cycle which 
necessitates uranium-235 (U235) enrichments as high as 93%. 
Obviously, to meet nonproliferation objectives, alternatives 
to this fuel cycle must be evaluated. 

The HTGR coated-particle fuel design allows a variety of 
alternative cycles, each having characteristic advantages and 
disadvantages in particular circumstances. Designing alterna- 
tive fuel cycles that can use the same core layout, fuel zoning 
and fuel blocks as the HEU/Th cycle means that conversion from 
one fuel cycle to another can be accomplished in available 
plant designs merely by replacing fuel blocks in each segment 
with alternative cycle blocks during the course of the regular 
refueling sequence. Since the parameters of f issile-to-f ertile 
and heavy metal- to-modera tor atom ratios control the conversion 
ratio of the reactor system, the versatility of the HTGR makes 
it a logical choice for tailoring fuel loadings and burnups, 
based on changing economic and resource availability conditions. 

Even more versatile and convenient for implementation of 
alternative fuel cycles among the advanced HTGRs is the pebble 
bed reactor [L4,Tl]. The injection-molded spherical graphite 
fuel balls used to fuel this system can easily be fabricated 
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with any desired fissile or fertile fuel density, and immedi- 
ately be introduced into the continuous refueling cycle. 

Operating an HTGR as a near- to-net-breeder utilizing the 
thorium fuel cycle requires low power density operation with 
heavy thorium loadings. Since operating trends are dictated 
by short term energy costs, obtaining a maximum conversion 
ratio and better fuel utilization is not necessarily a desirable 
near term goal. In obtaining a higher conversion ratio, in- 
creasing the thorium loading in the core while decreasing the 
core power density and fuel residence time would incur off- 
setting economic penalties which in turn may or may not be 
balanced by the advantages of improved fuel utilization. The 
penalties would depend on projected economic conditions, the 
projected availability of uranium ore, costs for separative 
work, reprocessing and recycle, and transportation and disposal 
costs . 

With increasing thorium loading, although an increase in 
the conversion ratio will be realized, a higher fissile inven- 
tory is also required. The resulting higher initial costs will 
be attractive if sufficient money is saved over the life of 
the reactor. This outcome is favored in an economic environment 
where large price increases are projected for fresh feed fuel. 

Lowering the core power density in an HTGR will result in 
a larger core volume for a given thermal output. Hence, addi- 
tional volume is available for fertile thorium loading, pro- 
ducing the same net result as previously described. Addition- 
ally, neutron utilization is improved due to a reduction in 
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losses to xenon-135 (Xel35) and protactinium-233 (Pa233) [K2]. 

The larger reactor volume will increase initial capital invest- 
ment and will also change the required coolant pumping power 
of the system. Keeping the coolant temperature increase con- 
stant across the core while maintaining similar height-to- 
diameter ratios, the larger, low power density core will require 
a higher coolant mass flow for a given thermal output. Because 
the HTGR uses helium as its primary cooling medium, very high 
coolant mass flow rates require that a significant fraction of 
the plant power output be utilized in operating the coolant 
pumps. Therefore, a variation in the core power density that 
will increase the pumping requirements can add a sizeable 
financial burden to the overall operating costs. 

In addition to the HEU/Th fuel cycle, the use of a uranium- 
233/thorium (U233/Th) cycle, a low-enriched uranium (LEU) cycle 
or various plutonium/uranium (Pu/U) or plutonium/uranium/ thorium 
(Pu/U/Th) fuel cycles are also possible in the KTGR [Bl,S3,Tl]. 
For all of the options the reduced total uranium ore requirement 
achieved by increasing the conversion ratio requires accepting 
either a higher capital investment or higher operating costs 
in other areas . 

The present work focuses primarily on the U233/Th cycle 
(because of its well known superiority in thermal and epi- 
thermal systems) and the possibility of its utilization in a 
system designed to accomplish most key nonproliferation 
objectives. Since U233 is required as the initial fissile 
loading in this cycle, the resource conservation is not direct. 
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but rather indirect. To provide the initial U233 loading, 
adequate fissile production in prebreeders, operating expressly 
for the production of U233 is required. To accomplish the 
objectives of nonproliferation effectiveness we will evaluate 
the feasibility of using the U233/Th cycle in either the 
internally recycled or the "Once-Through-Then-Out" (OTTO) 
modes in a pebble bed reactor system [Tl] . 

In the projection of any long-term energy strategy, the 
extension of uranium resources must be considered as most 
important, even in view of the increasing relevance of the 
objectives of nonproliferation. The present work, however, 
will attempt only to evaluate the feasibility of a single 
scheme for effectively providing a viable solution to the 
proliferation problem of safeguarding sensitive nuclear 
reactor materials through all steps of the fuel cycle without 
regard to its long-term effects on either the utilization of 
resources or the economics of implementing such a scheme 
commercially. However, because breakeven breeding is neces- 
sary to fully achieve the nonproliferation objectives, long- 
term ore utilization should be quite efficient. 

1 . 3 Previous Work 

The use of the breed/burn concept was proposed by Fischer 
et al [FI] for implementation in a reactor operating with a 
hard spectrum core. Depleted uranium monocarbide feed fuel would 
be introduced to breed fissile plutonium, and after a period 
of residence in the hard spectrum the fuel, now slightly 
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enriched, would be moved to the more thermalized outer region of 
the fast core. Energy would be generated as the fuel burns 
to a design-dependent burnup limit while the fissile- to-fertile 
fuel ratio is reduced toward the equilibrium value for that 
neutron energy spectrum. 

The present work extends the concept of the breed/burn 
cycle to assess the feasibility of its implementation in the 
thermal/epi thermal spectrum of the HTGR. Nonproliferation 
objectives of reducing fissile enrichments can be realized 
in addition to the benefits of decoupling this type of reactor 
from breeders and the constraints related to reprocessing 
facilities. The exception, of course, is the need to obtain 
the initial startup core, here requiring a substantial U233 
inventory . 

The introduction of this fueling scheme into the pebble 
bed reactor system provides the ability to readily adjust the 
core neutron spectrum simply by varying the modera tor-to— 
fuel (carbon-to-heavy metal) ratio. Additionally, core 
zoning (establishing separate breed and burn regions) can 
easily be accomplished in this application. A given ball, 
during its lifetime, would reside in two distinct regions. 

Radial zoning would make coolant temperature control difficult. 
Since the breed zone would operate at a lower power density 
than the burn zone, cold coolant bypass would degrade the bulk 
core exit temperature. Therefore, the core will be divided 
into axial zones, allowing it to be seen that for a fertile load 
ball, the "breed zone" can be defined as that region of the 



reactor in which the fissile concentration of the ball is in- 
creasing. Similarly/ the "burn zone" is that region in which 
the fissile concentration decreases. The desired breeding and 
burning effects can be accomplished by operating the breed zone 
in a predominantly epithermal spectrum. Thus, in its lifetime, 
a f ertile-f ueled ball would travel through both zones, enter- 
ing the thermal burn zone (together with pure graphite balls 
necessary to attain the desired moderator-to-f uel ratio) 
after first passing through the breed zone (which contains 
only fueled balls) . A schematic diagram of system operation 
is illustrated in Figure [1.1]. 

In successive passes through the breed zone, the concentra- 
tion of U233 is built up to the maximum possible amount for 
the given spectrum. It is then cycled (with additional gra- 
phite-only balls) into the burn zone, again being recycled 
for successive passes until a design-limited burnup is 
attained. At this point the depleted balls are removed from 
the system for disposal. 

Operating on this type of fuel cycle, the proposed reactor 
would allow very high burnups, and as a result, the total 
fissile enrichment of the discharged balls would be very low, 
particularly if the number of passes could be minimized, 
in which case discharge compositions approaching those of the 
OTTO cycle could be envisioned. If the reactor core can be 
adequately segregated so as to maintain distinct breed and 
burn regions, then the practicality of achieving the proper- 
ties sought by Fischer et al can be evaluated. 



Makeup Th232 feed 
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FIGURE 1.1: PEBBLE BED HTGR BREED/BURN SYSTEM SCHEMATIC 
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Work on HTGRs has been devoted mainly to the use of the 
HEU/Th fuel cycle, both in traditional graphite block-element 
reactors and pebble bed systems. Many alternative cycles, 
previously noted, have also been investigated more recently 
and in less detail, both from their ability to meet nonpro- 
liferation objectives as well as their respective roles in the 
extension of uranium resources. Current studies [D2,T1,V1] 
are focusing on the flexibility of the HTGR fuel cycle, which 
permits it to adapt to variations in fuel materials relatively 
efficiently. While favorable in many respects, this feature 
also makes it difficult to determine optimal strategies. This 
extra flexibility, relative to other reactor types, results 
from the particle loading system, which allows one to vary 
the fertile and fissile loadings almost independently and also 
to adjust the resonance shielding, not only by varying the 
size and loading of the fuel elements, but also by adjusting 
the dimensions of the individual fuel kernels [Dl] . 

Utilizing this flexibility it is possible to operate 
HTGRs with conversion ratios slightly greater than unity 
[R2,S3]. But considering resources and economics, the high 
fissile inventories required and the resulting recycle losses 
have made fuel designs with conversion ratios on the order 
of 0.9 appear more promising [S3], because near this value, 
burnup and processing losses reach a minimum [Dl] . 

Fuel cycle optimization in the HTGR has traditionally been 
most heavily concerned with minimizing fuel cycle cost. It is 
clear that this approach does not fully utilize the potential 
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in the HTGR for resource conservation through the use of the 
U233/Th fuel cycle. To obtain a breeding system with this 
cycle, if less than g-2 neutrons are lost to parasitic ab- 
sorbers, then after one neutron is utilized in sustaining 
criticality, one additional neutron is available for absorp- 
tion by the fertile thorium to produce, after B-decay, a 
new fissile U233 atom. The selection of the U233/Th cycle 
for this application is clear when the spectrum-averaged 
neutron yield per absorption, shown in Table [1.1], is con- 
sidered for this reactor type. For various carbon-to-heavy metal 
atom ratios (N^/N^) only fissile U233 has a value of n ade- 
quately high to allow breeding. 

The high flexibility in the design of fuel elements and 
the reactor refueling scheme in the pebble bed HTGR has made 
its use as a thermal breeder possible [R2] . Using highly 
enriched U233 ( (U233/U) - 0 . 9) produced in a prebreeder as feed 
fuel, thick radial blankets of thorium, high heavy metal load- 
ings in the fuel elements and relatively low heavy metal burnups 
will allow the pebble bed system to be self-supporting; inde- 
pendent of uranium resources after the first core inventories 
have been built up. Such designs have been shown to achieve 
conversion ratios greater than 1.05 [R2] . 

Other studies [T1,V1] , dealing more with nonproliferation 
and ore conservation strategies than with breeding potential 
have shown that proliferation resistance is best achieved in 
the low-enriched mixed oxide cycle (mixed uranium/ thorium 
oxides with 0235/0^0.20). Reactor loading is restricted to 
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TABLE 1.1 

SPECTRUM-AVERAGED NEUTRON YIELD 
PER FISSILE ABSORPTION* 

CARBON-TO-HEAVY METAL 
FISSILE ATOM RATIO (N C /N HM ) + 



ISOTOPE 


110 


250 


325 


450 


U233 


2.209 


2.240 


2.246 


2.253 


U235 


1.928 


1.989 


2.001 


2.011 


Pu239 


1.775 


1.788 


1.791 


1.794 



C/HM ratios correspond to spectrum averaged cross- 
section data given in Appendices B and D. 



t 
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non-sensitive materials while 95% of the bred fissile plutonium 
is burned in situ [Tl] . 

Resource conservation, however, is most favorable for OTTO 
cycles using highly-enriched U/Th mixed oxide ( (U235/U) -0.93) 
fuels designed for burnups of 100 to 130 GWD/MTHM [Tl] . 

From the above discussion, it can be seen that the wide 
latitude of compositions and fueling schemes available in the 
HTGR allow considerable breadth in designing strategies empha- 
sizing objectives of proliferation resistance, fissile breeding 
potential, uranium resource conservation or fuel cycle econom- 
ics. Each strategy has characteristic advantages and disadvan- 
tages, depending on the design criteria for the chosen fuel 
cycle. The breed/burn concept, when applied to the pebble bed 
HTGR system, combines the objectives of nonproliferation and 
fissile breeding. It is anticipated that the HTGR system can 
best provide the versatility needed to devise such a scheme for 
a thermal/epithermal system. 

1 . 4 Outline of Present Work 

The main body of this report consists of three chapters. 
Chapter 2 discusses the operation of the breed/burn fuel cycle 
in the pebble bed system and the relationship of the neutronic 
parameters of the two regions. Additionally, an evaluation is 
made of the potential of the pebble bed system for successful 
operation utilizing the breed/burn concept using one group 
models. Discussion is presented of the effects of spectrum 
hardening on the rate and amount of fissile breeding obtainable. 
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Chapter 3 describes the computational methods used for 
evaluation of the reactor system, nuclear data input, cross- 
section generation and the reactor model. 

The results of the analysis are discussed in Chapter 4 
with regard to reactor criticality, burnup limits and rates 
of fissile production as a function of composition. Possible 
modifications to the system are suggested. 

Chapter 5 summarizes the analysis and reiterates the main 
conclusions. Recommendations are made for future work. Appen- 
dices containing supplementary information are also included. 
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II. BREED/BURN OPERATION 



2 . 1 Introduction 

The pebble bed HTGR is expected to provide the best avail- 
able system for implementation of the breed/burn concept into 
a thermal/epithermal reactor because of it compositional 
flexibility and high burnup. In order to utilize the breed/ 
burn fuel cycle in this system, the feasibility of operating it 
as a two-zone reactor will be investigate. Dividing the 
core axially will allow the two zones to exist in harmony with 
the flow-through fuel cycle, such that all fuel elements will 
move through both zones during successive passes through the 
core. 

The purpose of this chapter is to consider the effects of 
various parameters on the rate of U233 production in the breed 
zone of such a system and to evaluate the potential for success- 
ful operation in the pebble bed HTGR. Possible modifications 
will be discussed to enhance system performance. 

2 . 2 Breed/Burn Potential for U233 Production 

The objective is to determine if, in fact, it is a feasible 
consideration to attempt using the breed/burn system in a peb- 
ble bed HTGR. For this analysis, the system considered will 
alternate breed and burn regions axially for the length 
of the reactor. The pebble bed HTGR is a very large, low 
power density system, and hence the net neutron leakage from 
the core is relatively low. By varying the composition 
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and geometry of each region, the respective neutrohic charac- 
teristics of the zones can oe altered, thereby changing the 
internal leakage between zones and hence changing the operating 
envelope of the entire system. 

First, an analytical method will be used to evaluate the 
parameters having the greatest effect on the rate of fissile 
breeding, and then the entire system will be evaluated for its 
ooerational capabilities. 



2 . 3 Factors Affecting the Breeding Ratio 

The customary index of fissile breeding (here of U233 
from neutron capture in fertile Th232) is represented by the 
breeding ratio, b, defined as: 



Rate of fissile production in the system 
Rate of fissile consumption in the system 



or. 



b = 



C 02 + C. 02 
c bz 

23 23 

A + A, 
c bz 



( 2 . 1 ) 



where C = neutron capture rate in the species indicated, 

A = neutron absorption rate in the species 
indicated, 

c = burn zone 

and bz = breed zone. 



For simplicity in this analysis, the effects of only the 



most significant materials will be considered. While it is 
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realized that a rigorous approach which would include the pro- 
duction and depletion of protactinium and U235 and the deple- 
tion of Th232 could yield more precise conclusions, the 
factors controlling the rate of fissile U233 production can 
be examined adequately by this more simple analysis. 

In each region a neutron balance can be written as: 



Production = Losses + Absorption + Leakage 



which, for the burn zone, is given by: 



VF 23 + VF 02 = L + F 23 + F ° 2 + C 23 + C ° 2 + A P . (2.2) 

c C cccccc 



Similarly for the breed zone: 



vF bz 23 + vF bz 02 + L c “ F bz 23 + F bz 02 + C bz 23 + C bz 02 + A bz P (2 - 3) 



where 

F = fission rate in the species indicated, 

L = neutron leakage from the region indicated 

P 

and A = parasitic absorption. 



Combining these terms the breeding ratio of the system can be 
expressed as: 



(v - 1) F, 



02 



A 



“ n s + 



A 



23 



23 



b 



1 



(2.4) 
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where n s 



mean neutron yield per fissile absorption in the 
system 



v 




A 23 
s 



(2.5) 



and s sum of indicated interactions in the breed 

and burn zones . 



Similarly, the breeding ratios for the burn zone alone 
and the breed zone alone, respectively, are given by: 



b . n + lv -» F c° 2 . 

c c A 23 ' A 23 - 1 



and 



(v - 1 F 02 A P 

b « n + _ , 

bz 'bz 23 , 23 X 



where 



— 23 

V F J 

= 9 

c - 23 

A 

s 



v F 



23 



and 



n 



bz 



bz , 23 

A s 



( 2 . 6 ) 



(2.7) 



( 2 . 8 ) 



(2.9) 



Equations (2.4), (2.6) and (2.7) are composed of three main 

terms. The first term, g , dominates the expressions, and 
hence increasing g will increase the breeding ratio proportion- 
ately. This can be accomplished by hardening the core neutron 
spectrum through an increase in the concentration of the 
fertile and fissile species. That is, lowering the carbon-to- 
heavy metal (C/HM) atom ratio and thus decreasing the modera- 
tion results in a harder neutron spectrum. 
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The second term, (v - 1) F^/A^, represents the gain in 
D due to fast fissioning of fertile material. Because of the 
small threshold fission cross-section of Th232 this term is 
relatively small. 

P 23 

The third term, A /A , represents the neutrons lost to 
parasitic absorption and leakage and hence is negative. Thus, 
the objective of increasing the amount of fissile breeding 
necessitates minimizing the amount of parasitic absorption 
and leakage losses. Fuel material absorption cross-sections 
and those of structural materials, therefore, become important 
design considerations. From this analysis we must conclude 
that in order to increase the breeding ratio in our system a 
hardened neutron spectrum is generally favored since it mini- 
mizes parasitic absorption and maximizes neutron yield, as will 
be shown in succeeding sections. Leakage losses may or may 
not be improved, however. The first design strategy clue is 
then to obtain as hard a spectrum as practical in both the 
breed and burn zones within other constraints. Because we 
want to be able to produce adequate masses of fissile material 
in the breed zone to fuel the burn zone, our attention 
focuses on the parameters of the breed zone. It is realized 
at this point, however, that the two zones will have significant 
coupling in their characteristics. That is, considering a 
given number of neutrons in the system, favoring neutron avail- 
ability in one region will decrease the effectiveness of the 
other. If the breed zone is, under any conditions, unable to 
breed adequate masses of U233, then the neutron balance in the 
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burn zone will still be unable to effect satisfactory system 
operation. It is for this reason that the breed zone para- 
meters will be examined in more detail. 

2 . 4 Breed Zone U233 Conversion 

To some degree, breeding of U233 will occur throughout 
the reactor. In applying the breed/burn concept to the pebble 
bed HTGR a criterion for successful operation has been set 
that adequate breeding in the breed zone alone can be achieved 
to provide sufficient total fissile mass to fuel the burn zone. 
Obviously, a simultaneous objective is to maintain system 
criticality for a reasonable cycle lifetime. Breeding of U233 
in the entire system can be divided into two parts, correspond- 
ing to burn zone breeding and breeding in the breed zone. 

For the burn zone: 



b = 



Rate of fissile production in the burn zone 

Rate of fissile consumption in the breed and burn zones 



02 



(A 23 + A h 2 V 
c bz 



( 2 . 10 ) 



while for the breed zone: 



_ Rate of fissile production in the breed zone 

) bz ” Rate of fissile consumption in the breed and burn zones 



02 



'bz 



23 

(A + A, 
v c bz 



23 , 



( 2 . 11 ) 
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Solving Eq. (2.3) for the neutron capture rate in Th233 
in the breed zone: 



C 02 = 7 F„ 23 + VF. 02 — F, 23 -F k 02 - Cv , 23 - A. P +L 
bz bz bz bz bz bz bz c 



= (v - 1) (F 23 + F ° 2 ) - C. 23 - A P + L . 

bz bz bz bz c 



Using 



F = F 23 + F 02 and A, PC = C, 23 + A, P , 
bz bz bz bz bz bz 



we can write: 



(v-l)F, + L -A, 

_ bz c bz 

bZ_ <A c 23 ^ bz 23 > 



PC 



( 2 . 12 ) 



Since the thermal power of an operating reactor is essen- 

23 23 

tially constant, then A + A. , the total rate of fissile 
u c bz 

consumption in the reactor, can be treated as a fixed value. 

It can be seen then, that Eq . (2.12) expresses the dependence 

of the rate of U233 conversion in the breed zone on: 

(a) the rate of neutron leakage from the burn region 
into the breed zone, [F> c ] , 

(b) the rate of neutron production in the breed zone, 

[ (v - 1) F bz ] and 

PC 

(c) the rate of neutron loss in the breed zone, [A bz ] • 

Losses due to parasitic absorption and leakage in the 
breed zone have the effect of cancelling the gains made by 
the production of neutrons in the zone. It is these factors 
that must be reduced while simultaneously increasing the rate 
of neutron production in the breed zone to secure high neutron 
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availability and improve the rate of U233 conversion, 
Combining Eqs. (2.10) and (2.11) yields: 



C 02 + C 02 
b c + b bz = c C bz 

c bz 



(2.13) 



which, of course, is the system breeding ratio, expressed in 
Eq. (2.1). It is clear then, that the rate of breeding in the 
breed zone is increased only at the expense of breeding in the 
burn zone. The objective then becomes to optimize the total 
rate of U233 conversion for the system, under the constraints 
that criticality be maintained -- which will be made more diffi- 
cult if the burn zone breeding ratio is small. Since this breed/ 
burn concept is designed for application in commercial power 
plants, high burnups and good ore utilization are essential. 
Design parameters that will allow the production of U233 in the 
required quantities while achieving reasonable fuel element 
lifetimes and core power densities must, therefore be considered. 



2 . 5 Evaluation of Homogeneous Zone Characteristics 

Operation of the breed/burn fuel cycle for the pebble bed 
reactor system has been defined and its operation described in the 
preceding sections. To predict more precisely the characteristics 
of such a system, the parameters involved in breed/burn operation 
for homogeneous one-zone regions will be evaluated. This analysis 
should show, for a given composition, how a region would operate 
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without the coupling effects (leakage) of a zone of different 
composition on its boundaries. 

Neutron economy in any reactor system is a function of 
the ambient neutron spectrum, axial and radial zoning of 
compositions, power level and system geometry. The rate of 
neutron leakage, being dependent on the size of the reactor, 
decreases with increasing core volume, holding other para- 
meters constant. For this analysis, to provide the most 
ideal and hence a limiting-case representation of the feasi- 
bility of breed/burn operation, an infinite-medium model will 

be used. Pebble bed reactors, designed in cylindrical geo- 

3 

metries, operate at low power densities (typically 5 w/cm 
[W3] ) and hence are inherently large systems. Reactors having 
thermal power outputs currently competitive commercially 
(v3000 MW would therefore have core volumes on the order of 

3 

600 m , and in such cases would have neutron leakage fractions 
below 5% [W3] . Using an infinite system model for this analy- 

sis, therefore, is not an unreasonable approximation. 

For layered zoning of the pebble bed system described in 
Chapter 1, both breed and burn zones will be treated here as 
uniform, separate regions. In this way, for a given spectrum, 
the best possible operating characteristics can be determined 
as there will be no coupling effects between regions which 
would perturb the spectrum in the vicinity of the boundary . 

As shown in the previous section, the breeding ratio for the 
system is simply the sum of the breeding ratios in the breed 
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and burn zones. Therefore, if one can determine the maximum 
amount of breeding possible in a homogenized region, this will 
approximate the system-averaged breeding ratio for a two-zone 
reactor in the ideal case. 

The pebble bed HTGR utilizes injection-molded spherical 
graphite elements, illustrated in Fig. [2.1], fueled with 
pyrolytic carbon (PyC) and silicon carbide (SiC) coated ThO^ 
and UC >2 microspheres. The balls themselves provide all of 
the fuel, moderator, structural material and coolant channels 
(free space between balls) necessary for the reactor. There- 
fore, in a clean core, the only materials present are those 
found in the newly fabricated feed balls, in addition to the 
helium coolant. As a function of the carbon- to-heavy metal 
atom ratio, the number densities of fertile feed ball materials 
are given in Table [2.1] (for calculation of number densities, 
refer to Appendix D) . 

Variation in the carbon- to-heavy metal atom ratio signifi- 
cantly alters the neutron spectrum in a given region due to the 
change in moderation. One-group constants for carbon-to-heavy 
metal atom ratios of 110, 250, 325 and 450 are given in Appen- 
dix A, collapsed from the AMPX-generated 4-group cross-sections 
given in Appendix C (the AMPX system is described in Chapter 3) . 
To illustrate the effect of the carbon-to-heavy metal atom 

ratio on the neutron spectrum, the spectrum- averaged values 

23 

for n , the rate of neutron production per fissile absorption. 



given by: 
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FIGURE 2 . 1 



PEBBLE BED HTGR FUEL ELEMENT 
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TABLE 2.1 



COMPOSITION* OF PEBBLE BED FEED BALLS 



C/HM ATOM RATIO 



NUCLIDE 


110 


250 


325 


450 


Thorium 


2.422E-4 


1 . 04 lE-4 


7.872E-5 


5.702E-5 


Oxygen 


4.844E-4 


2.083E-4 


1.574E-4 


1.140E-4 


Carbon 


2 . 666E-2 


2. 610E-2 


2. 600E-2 


2.591E-2 


Silicon 


6 . 644E-4 


2.857E-4 


2. 160E-4 


1.564E-4 


* Number 


densities 


-1 -1 
in b cm 
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vcr 



23 



23 ' 
a 

a 



(2.14) 



the microscopic fission cross-section for Th232, 



02 



and the 



microscopic neutron absorption cross-sections for the parasitic 

absorbers, a carbon , a oxygen and a silicon^ are tabulated in 

d a. 

Table [2.2] . 



It can be seen from these results that the preference for a 

hardened neutron spectrum is not immediately justified. The 

values given in Table [2.2] are those used in determining the 

breeding ratio of the system, expressed by Eg. (2.1). It is 

immediately apparent that increasing the carbon- to-heavy metal 

atom ratio from 110 to 450 increases n . This is opposed to 

23 

our analysis m that an increase m n with added moderation 

would cause an increase in the breeding ratio. However, Table 

[2.2] also reveals decreases with a harder spectrum in the 

absorption cross-sections of the parasitic nuclides and an 

increase in the fission cross-section of Th232. These factors 

lead to increases in the breeding ratio with spectrum harden- 

2 3 

ing. In fact, while ri decreases a- 2 % from C/HM = 450 to 
C/HM = 110, the macroscopic parasitic absorption cross-section 
(Z + E ^ + Z ^ 1 ) decreases ^36%, while the Th232 fission 

cl cl cl 

cross-section increases ^72%. So although Eq. (2.1) is domi- 

23 ... 

nated by n > the much larger changes in Th232 fissioning 

and parasitic absorption could have the dominant effect on 

the breeding ratio. For these reasons it is important to 

more carefully consider the effects of the ambient spectrum 

on breed/burn operation. 
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TABLE 2 . 2 

SPECTRUM-AVERAGED HOMOGENEOUS ZONE PARAMETERS 



C/HM ATOM RATIO 



PARAMETER 


110 


250 


325 


450 


n 23 


2 . 209 


2 . 240 


2.246 


2.253 




8 . 097E-3 


6. 051E-3 


5. 446E-3 


4.718E-3 


a C 

a 


3. 830E-4 


7 . 095E-4 


8.191E-4 


9. 507E-4 


Q 

S» O 


5. 306E-4 


4 . 146E-4 


3.811E-4 


3 . 409E-4 


a Si 

a 


1.358E-2 


2. 709E-2 


3.160E-2 


3 . 702E-2 



Cross-sections in barns. 
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2 . 6 Variation in Neutron Spectra 

The proposal by Fischer et al [FI] for implementation of 
the breed/burn fuel cycle in a hard spectrum core indicated 
that, with respect to breeding ratio, the best results would 
be obtained if the ambient spectrum was hardened as much as 
possible. In considering the pebble bed HTGR, our first indi- 
cations verify Fischer's claims for the preference of the 
hardened spectrum. It has now become necessary to justify 
this claim more fully if it is to be used as a design objec- 
tive for thermal reactor versions of the breed/burn concept. 

Two criteria for successful operation of the breed/burn 
cycle have been established: (1) the ability to maintain 

system criticality, and (2) the ability to breed adequate 
fissile mass in the breed zone to fuel the burn zone. The 
second objective is solely a function of the breeding ratio in 
the breed zone. Maximizing the rate of U233 conversion by 
increasing the parameters in Eq . (2.1), as described in the 

previous section, makes it easier to achieve the second goal 
alone, if adequate exposure time can be provided by artifi- 
cially satisfying the first objective. However, when the 
criteria for successful operation must be satisfied simul- 
taneously, it is realized that the optimization of one 
parameter detracts from the ability to satisfy the other. 

By this reasoning, Fischer et al, in achieving both 
breed/burn operation and the highest breeding ratio possible, 
might not be able to maintain criticality in their overall 
system. If this is the case, operating a breed/burn cycle 
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becomes a futile task. Then, not only must the factors in- 
volved in obtaining a high breeding ratio be considered, but 
also those factors affecting the criticality of the system. 

Sheaffer et al [SI] have correlated microscopic cross- 
sections with the ambient neutron spectrum in fast reactors 
by defining a spectral index parameter, S, given by: 



S 



v£ 



f tr 



(2.15) 



Physically, S is a measure of the average neutron energy in 
the core population. Cross-sections, as a function of S, are 
given by: 



a = 




(2.16) 



where and g are constants appropriate for the indicated 
nuclide and cross-section (the constants and g are tabula- 
ted in Appendix B, Table [B.l] for the materials used in this 
analysis). The correlation expressed by Eq. (2.16) was devel- 
oped for application to fast reactor systems and has been 
verified (with very good accuracy) for values of S in the 
range 0.3 to 1.0 [SI] (corresponding to a wide range of fast 
reactor spectra, from oxide-fueled LMFBRs to metal-fueled 
GCFRs) . Sheaffer's method can be used to evaluate breed/burn 
performance in FBR systems . 
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The ability of the system to breed fissile enrichments 
suitable to fuel its fast core must be investigated. In any 
breeding system, if it was assumed that the neutron flux could 
be maintained artificially, then the system would eventually 
achieve a conversion ratio of 1.00, at which time the rate 
of fertile neutron capture would exactly equal the rate of 
fissile destruction by neutron absorption. Thus, the system 
would maintain an asymptotic fissile enrichment limit which 
would be the maximum possible enrichment compatible with the 
ambient spectrum. 

If the fertile capture rate is equal to the fissile ab- 
sorption rate: 



No cf> 
c 



fertile 



No <j> 
a v 



fissile 



(2.17) 



Neglecting all heavy metals except the principal fertile and 
fissile species, the fissile enrichment, e, can be defined as: 



„f issile 
N 



„,fissile , ..fertile 
N + N 



(2.18) 



Combining Eqs . (2.17) and (2.18), an expression for the 

asymptotic fissile enrichment, £ AS y/ as a function of the 

, . fertile, fissile . 
cross-section ratio a /a is: 

C 0. 



'ASY 



fertile, fissile 
a /o 

c a 

n , i fertile, fissile 
l+o /0 = 



) ■ 



(2.19) 
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For a plutonium- 2 3 9/uranium- 23 8 system, Eq . (2.19) becomes: 



„ 28, 49 

o /o 
c ' a 



"ASY , , 28, 49, 

1 + <°c /0 a 1 



and similarly for a U233/Th system: 



„ 02, 23 

o /o 
c ' a 



'ASY , , . 02, 23, 
1 + (a /a ) 
c z 



( 2 . 20 ) 



( 2 . 21 ) 



With these approximations the practical end-of-life (EOL) 
fissile enrichment, e , obtainable with high fuel burnup 
in a breed/burn system, can be roughly estimated. 

Consider the following conceptual picture of breed/burn 
operation: pairs of physically identical fuel elements exist 

— one in the breed zone and the other in the burn zone. The 
fissile enrichment in the breed zone at end-of-cycle must 
equal the fissile enrichment in the burn zone at the beginning- 
of-cycle if a "stand-alone" breed/burn cycle is to be feasible. 
For every fission in the burn zone element roughly one excess 
neutron is produced which can be used to produce fissile 
material in the breed zone element. Since both elements have 
very nearly equal heavy metal contents, the maximum percent 
fissile enrichment in the breed element is equal to the percent 
burnup in the burn zone. That is, for breed/burn operation 

(in %) : 

B c [MWD/MTHM] 

£ BOL “ e EOL, ~ 774 

C bz 10 



( 2 . 22 ) 
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The potential of different reactor types is as follows: 
LWR: e BOL^ 3%; B EOL ,V 30,000 MWD/MTHM =3% 

HTGR: e BOL a ' 9%; ^OL^ 120,000 MWD/MTHM = 12% 

LMFBR or GCFR : e BQL ^15%; B EOL ^ 100 ' 000 MWD/MTHM = 10%. 

This information is illustrated in Fig. [2.2], where it can 
be seen that the potentially feasible breed/burn systems would 
be those operating below the line, whereas those above the line 
would not be feasible. The circumstances indicated in Fig. 

[2.2] are, in fact, a large part of the motivation for investi- 
gating the HTGR as a breed/burn reactor in the present study. 

Since the cross-sections of the fertile and fissile 
species vary with the ambient spectrum, Eq . (2.19) can be 

used to estimate the attainable asymptotic enrichment limit 
for any given system. This limit indicates for a reactor 
operating with the breed/burn fuel cycle, whether to harden or 
soften the ambient spectrum. A high asymptotic enrichment 
would imply that relatively short burnups would be required 
to achieve a desired enrichment value, whereas a low asymptotic 
limit would mean that if a beginning-of-lif e enrichment was 
required near the asymptotic value, then very high burnups 
would be essential to achieve that limit. This argument adds 
preference to those spectra whose asymptotic enrichment limits 
are as high as possible. 

From Eq. (2.16), cross-sections can be determined for any 
given fast reactor spectrum (composition). For this analysis. 
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FIGURE 2.2: SCHEMATIC DIAGRAM SHOWING POTENTIAL 

FEASIBILITY OF BREED/BURN SYSTEMS 



*Typical FBR with current burnup limits. 
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the following approximation is used: 



a fertile = fertile (2.23) 

c a 



which, from Eq. (2.19), yields: 



'ASY 



fertile, fissile 

q a /° a 

1 + (a fertiie/o fissile, 
3. a. 



(2.24) 



This approximation overestimates the value of e ASy and there- 
fore yields a conservative upper limit for the asymptotic 
enrichment. Listed in Table [2.3] are values .for e A g Y ' us i n 9 
Eq . (2.24), calculated for spectra varying from S = 0 . 3 to 

S=1.0. Also tabulated are values for fissile n, given by: 



n 




(2.25) 



for the indicated fissile species. It is clear from the data 
in Table [2.3] that values for fissile n for both the Pu239/U 
and U233/Th systems increase with spectrum hardening, justify- 
ing the claim of Fischer et al that harder spectra will yield 
better breeding potential. This alone, however, does not 
justify the preference of the harder spectra for the breed/ 
burn concept. In fact, because a balance must exist between 
fissile breeding and system criticality, the asymptotic fissile 
enrichments, illustrated in Fig. [2.3], indicate that the 
softer fast reactor spectra would be preferable. Operation 
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TABLE 2.3 

ASYMPTOTIC ENRICHMENTS AND FISSILE n VALUES 
FOR FAST SPECTRA* 



s 


£ 2 3 
ASY 


n 23 


£ 4 9 
ASY 


4 9 

0 


0.3 


0.1388 


2.175 


0.1480 


2.308 


0.4 


0.1210 


2.224 


0.1368 


2.415 


0.5 


0.1085 


2.264 


0.1286 


2.502 


0.6 


0.09918 


2.296 


0.1223 


2.574 


0.7 


0.09190 


2.324 


0.1172 


2 . 638 


0.8 


0.08591 


2.348 


0.1128 


2.694 


0.9 


0.08097 


2.370 


0.1091 


2.744 


1.0 


0.07675 


2.390 


0.1059 


2.790 



* From Sheaffer et al [Si] correlation: cross-sections 

are listed in Table [B.2] . 



ASYMPTOTIC ENRICHMENT 
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FIGURE 2.3: ASYMPTOTIC FISSILE ENRICHMENT 

LIMITS FOR U233 AND Pu239 SPECTRA 
USING SHEAFFER'S CORRELATION [SI] 
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can, in this case, achieve much higher fissile enrichments for 
a given burnup than can a system operating in a much harder 
spectrum. For a Pu239/U fueled system, the softer (S = 0.3) 
spectrum burnup (assuming a linear relation between burn zone 
and breed zone enrichments) would require an exposure time ■v 
32% less than the hard (S=1.0) spectrum, while for the U233/ 

Th system the same change in spectrum would reduce the exposure 
time by 'v 4 8%. 

These cases have indicated, contrary to the proposal of 
Fischer et al , that a softer spectrum is preferred for breed/ 
burn operation in a fast reactor. A similar analysis can be 
performed to determine the more favorable operating conditions 
for the thermal/epithermal HTGR . Because Sheaffer’s correla- 
tions were developed for fast reactor calculations they do not 
account for the high degree of resonance shielding in the 
fertile species for highly epithermal applications. Hence, 
the model is invalid for values of S< 0.3; the lowest values 
verified . 

A similar spectral index can be derived for application 
in the U233/Th HTGR cases. In the expression for the reso- 
nance escape probability from classical continuous slowing 
theory, p, given by: 



E dE‘ 
a 



e' e(£ a + i s > 



r 



p = exp 



(2.26) 
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an effective E^ can be used to account for leakage by adding 
a "leakage cross-section": 

E + E a + DB 2 . (2.27) 

d ci 

Then, accounting for criticality by applying the neutron 
balance: 



E^ + DB' 

ci 



= vE, 



(2.28) 



one gets: 

p = exp 



/ 



vE, 



K (vE f + E s ) 



dE' 

E ’ 



(2.29) 



Since E - E, 
s tr 

mal region vEf 
by : 

S ' 



, and recognizing that in the thermal / epither- 
« E^ r , a spectral index, S' can be defined 




(2.30) 



Using this definition, the index can be redefined as: 



S 



S ' 

1 + S' 



(2.31) 



which yields the same form used by Sheaffer et al [SI] , given 
by : 



S 



vE, 



vE f + ?E tr 



(2.15) 
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From the computations to be discussed in Chapter 3, values 
used to calculate S were found for the pebble bed HTGR composi- 
tions used in the analysis of Section [2.5]. Table [2.4] lists 
values for e ^gy / given by Eq. (2.21), with corresponding val- 
ues of S as a function of the carbon- to-heavy metal atom ratio. 
In Fig. [2.4] the asymptotic enrichments, ’ are Pi° tte< ^ 

versus the spectral index for both the U233/Th fast reactor 
cases and the pebble bed systems. It can be seen that in the 
thermal/epi thermal region, spectral hardening is advantageous 
in increasing the asymptotic enrichment limit. 

As shown previously, although the harder HTGR spectra 
lead to increasing values of fissile n» decreasing parasitic 
absorption rates and improved fertile fission rates combined 
with the benefit of higher asymptotic enrichment limits indi- 
cate a strong preference for the use of a hardened spectrum 
for successful thermal/epithermal breed/burn operation. This 
conclusion is contrary to that of fast reactors, where the 
softer spectra are preferred. The indication here is that 
an optimum breed zone ambient spectrum may be attainable in 
the regions between those considered in this analysis. The 
apparent maximum value for s ^ own i n Fig. [2.4], clearly 

indicates this possibility. Hence, in order to maximize the 
effectiveness of the breed/burn concept it would be advisable 
to consider systems operating with neutron spectra in the 
area of maximum e^gy between a "soft" FBR and a "hard" HTGR. 
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TABLE 2.4 

SPECTRAL INDEX AND ASYMPTOTIC ENRICHMENT 
LIMITS FOR THE U233/Th PEBBLE BED HTGR 



PARAMETER 


C/HM ATOM RATIO 
110 250 325 450 


S 


0.0724 0.0583 0.0518 0.0438 


£ asy 


0.03526 0.02686 0.02497 0.02302 



ASYMPTOTIC U233 ENRICHMENT 
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FIGURE 2.4: ASYMPTOTIC ENRICHMENTS FOR U233/Th 



SYSTEMS 
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2 . 7 System Criticality 

In the previous section it was shown that in using a 
breed/burn fuel cycle in the pebble bed HTGR, it is advanta- 
geous to harden the breed zone neutron spectrum in order to 
increase the rate of fissile U233 production. The required 
burn zone burnup for each element necessary to achieve any 
given breed zone fissile enrichment can therefore be shorten- 
ed, and hence, the number of passes through the breed zone 
required of each fuel ball can be minimized. 

The fissile enrichments attainable for each composition 

can be approximated using Eq. (2.21), again utilizing the 

one-group constants given in Appendix A. As a function of 

the carbon- to-heavy metal atom ratio, these enrichments are 

listed in Table [2.5] along with the cross-section ratio, 

02 23 

o /o . Using these values to approximate the beginning- 
of-life fissile enrichment for a clean core, the system criti- 
cality at the beginning of each cycle can be estimated. Apply- 
ing a one-group model, the infinite-medium multiplication 
factor, k^ , can be estimated at the beginning-of-lif e as: 



k 



f 



(2.31) 



00 



BOL 



Z 



a 



and the BOL excess reactivity, P B0L / is: 



k 



1 



oo 



BOL 




BOL 



(2.32) 
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TABLE 2 . 5 

ASYMPTOTIC FISSILE ENRICHMENT LIMITS 



PARAMETER 110 

23 0.03655 

a 

a 

e ASY (%) 3 ' 526 



C/HM ATOM RATIO 
250 325 

0.02760 0.02561 0 

2.686 



450 

02357 



2.497 



2.302 
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Values for and p are given in Table [2.6]. 

BOL BUij 

It can be seen that in all cases, k^ >1.00, indicating 

°°BOL 

criticality at the start of the cycle. Using the amount of 

U233 bred, for example, in the region of C/HM = 110 to fuel 

the regions where C/HM = 250 or higher would yield an even 

greater value for k^ . For the four compositions consid- 

BOL 

ered, the values of k TO obtained for various combinations of 

BOL 

breed and burn zone carbon- to-heavy metal atom ratios are 
shown in Table [2.7] . The highest k obtained is for that 

OG 

BOL 

of a burn zone with C/HM = 450 being fueled by a breed zone 
where C/HM = 110. This combination has the softest burn zone 
spectrum operating with the hardest spectrum (and hence, high- 
est e ASY ) breed zone of those compositions evaluated. 

Careful analysis is required in these cases to insure that 
a given breed zone would produce not only adequate fissile 
enrichment, but also sufficient total fissile mass to fuel the 
burn zone. From this analysis, the first order estimate leads 
to the conclusion that it is possible to breed enough fissile 
material in a breed zone, given adequate exposure time, to fuel 
and obtain criticality in a pebble bed breed/burn system. The 
amount of BOL excess reactivity is quite small, however, which 
indicates that it may not be possible to maintain criticality 
long enough to breed fissile enrichments near the asymptotic 
limits. One can adapt to this situation to some extent by 
using rapid fuel flow through the core, approaching continuous 
refueling in the limiting case, but of course this increases 
the complexity of operation. 
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TABLE 2 . 6 

BEGINNING-OF-LIFE INFINITE-MEDIUM 
MULTIPLICATION FACTORS 



C/HM ATOM RATIO (both zones) 



PARAMETER 


110 


250 


325 


450 




1. 435E-3 


7.782E-4 


6.163E-4 


4. 662E-4 


I 


1. 293E-3 


7. 110E-4 


5.695E-4 


4.392E-4 


a 


k 


1.110 


1.095 


1.082 


1.062 


°°BOL 


P BOL 


9 . 91 


8.64 


7.60 


5.78 
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TABLE 2 . 7 

FOR VARIOUS ZONE COMBINATIONS* 







BURN 


ZONE 


C/HM ATOM 


RATIO 






110 


250 


325 


450 




110 


1.110 


1.256 


1.286 


1.313 


BREED ZONE 


250 


0.9515 


1.095 


1.125 


1.152 


C/HM ATOM 


325 


0.9103 


1.052 


1.082 


1.109 


RATIO 


450 


0.8654 


1.005 


1.035 


1.062 



* in burn zone using asymptotic enrichments 
achievable in breed zone 



-60- 



The analysis of Section [2.6] can be refined further to 
consider some of the specific features of the breed/burn cycle. 
In the expression for , given by Eq . (2.31), a "leakage 

cross-section" can be added to the absorption cross-section, 
so that: 



vZf 

X a + S L 



(2.33) 



where, for example, one can use Z = DB Z . The criterion for 

Xj 

criticality, that k>^1.00, yields the necessity that: 



or. 



k > 1 + ^ 
oo — Z 



1 > 



23 



—23 02 P 

Z + Z + Z 
a a a 



(2.34) 



(2.35) 



where E refers to neutron absorption by the parasitic ab- 
sorbers . 

Therefore, Eq . (2.35) can be used to find that: 



k - 1 > n ? — T i p (2.36) 

1 + b„ + (Z ° 2 /Z 23 ) + (Z = P /Z_ 23 ) 
i a 



a a 



where b c is the breeding ratio in the burn zone, approximated 



by : 



02 



b c „ 23 



(2.37) 
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In Eq . (2.37), both £ ^and I P /Z ^ are small. The 

-La a a 



criticality limit becomes: 




(2.38) 



c 



In the burn zone, however, the maximum value of b is attained 

c 

under EOL asymptotic conditions. In this case, b c = 1.0, 
which therefore defines a lower limit for k of: 



oo 



k > 1.5 , 

OO ' 



(2.39) 



which is now a third criteria for successful system operation. 
Thus, k^ for the burn zone must exceed 1.5 for a breed/burn 
cycle to be able to operate independent of additional external 
fissile sources. 



Table [2.7], was k^ ^ 1.3, we must conclude that for the pebble 
bed HTGR, stand-alone breed/burn operation does not appear 
possible . 

2 . 8 Summary 

In this chapter the parameters affecting successful breed/ 
burn operation in the pebble bed HTGR have been examined. While 
preliminary evaluation suggested that spectral hardening would 
detract from system performance due to the decreasing values of 
fissile g , closer examination revealed that while ri decreases 
slightly, parasitic absorption decreases significantly while 



Since the greatest value obtained for k 



given m 



BOL 
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the fertile fission rate also increases substantially. In 
addition, harder spectra yield higher asymptotic enrichment 
values, thus making it easier to achieve acceptable fissile 
enrichments in the breed zone with reasonable burnups in the 
burn zone. The results also showed that highly thermalized 
systems are best for the burn zone because of their high k^ 
at a given load enrichment and low asymptotic (hence discharge) 
fissile enrichment. Of all reactors of current interest, 
it was found that the HTGR had the most favorable ratio of 
discharge burnup ro reload enrichment. 

These observations, with respect to spectral hardening, 
were found to apply for thermal and epithermal systems only. 

In fact, for fast reactors it was found that the breed function 
of breed/burn cycles would be best carried out in softer 
spectra, such as those found in oxide-fueled LMFBRs -- or, 
indeed, even softer. 

In examining system criticality, although it was found 
that a breed zone in a pebble bed reactor could fuel and 
achieve criticality (k = 1.0) in an infinite burn zone, the 
high BOL excess reactivity (k OT ^ 1.5) necessary for finite 
coupled zone (high burn and breed zone leakage) system opera- 
tion could apparently not be achieved. Hence, operating such 
a system in the pebble bed HTGR would require some degree of 
fissile makeup to operate for any sustained period; either 
that or a substantially harder breed zone spectrum (yielding 
higher asymptotic enrichment and hence larger k^ in the burn 
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zone) . An even higher carbon-to-heavy metal atom ratio in the 
burn zone could also be contemplated. 
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III. COMPUTATIONAL METHODS 

3 . 1 Introduction 

The analysis in the previous chapter led to the conclusion 
that in a pebble bed HTGR, it would probably not be possible 
to sustain a stand-alone breed/burn system. However, the one- 
group calculations performed could only indicate the gross 
features of system operation and were limited to either the 
beginning-of-lif e or the asymptotic limit of ultra-high burnup. 
The effects of fission product buildup, U235 and Pa233 produc- 
tion, Th232 depletion, fast fission, power density and group- 
wise changes in flux and cross-sections, all as functions of 
burnup, were neglected. When fully evaluated, the sum of the 
above effects can significantly alter the reactor performance 
characteristics as a function of time, which therefore neces- 
sitates a more precise evaluation. 

Calculation of the performance characteristics of the 
proposed breed/burn system, using state-of-the-art computer 
programs, is described in this chapter. The system model, 
cross-section generation, nuclear data input, codes used and 
problem descriptions are described. 

3 . 2 System Model 

In the breed/burn fuel cycle, it is necessary for the 
reactor system to operate under "breakeven" breeding condi- 
tions. That is, as a minimum, one fissile atom must be pro- 
duced in the breed zone for every one destroyed in the burn 



-65- 



zone. A first approximation can be used to evaluate this 
ability (i.e., by considering a homogeneous system) as was 
done in Chapter 2. The two- zone layered system, described in 
Chapter 1, will operate with a strong neutronic coupling be- 
tween zones. The burn zone fission rate, and hence, power 
density, will be higher than that in the breed zone. Also, 
the total neutron flux in the burn zone should be higher, 
whereas groupwise flux will be shifted towards the faster groups 
in the breed zone, assuming that the carbon-to-heavy metal atom 
ratio in the breed zone is less than that of the burn zone. 

A conservative overestimation of the performance of the 
system can be found, therefore, by homogenizing the two zones 
into a single region. If the conditions required for success- 
ful breed/burn operation can be met in the homogeneous region, 
then a more careful analysis is required to determine, when the 
reactor is modeled in a heterogeneous, axially-layered configura- 
tion, if satisfactory performance is still achieved. 

In analyzing the heterogeneous two-zone system, the re- 
quirement for the breeding of adequate total mass in the breed 
zone to fuel the burn zone dictates size limitations on the two 
zones, contingent upon the number of passes through each zone 
required to achieve design burnup limits. That is, given the 
compositions of each zone, in a period of one cycle, the total 
fissile mass being bred in the breed and burn zones must equal 
the total mass destroyed in the burn zone for net breeding; 
for breed/burn operation the breed zone alone must provide the 
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burn zone fissile deficit. Therefore, calculations must be 
performed iterating between size options for the two zones. 

If the reactor is assumed to be infinite radially, then the 
calculations can be performed using an infinite slab model 
varying only the compositions and axial thicknesses of each 
zone . 

The first analysis to be performed, however, utilizes a 
homogeneous zone as previously described. Unity albedo boun- 
dary conditions are used for all three dimensions, thus effec- 
ting an infinite reactor calculation. Again, an infinite slab 
treatment is appropriate. As will be described in Section [3.4], 
compositions (carbon-to-heavy metal atom ratios and fissile 
enrichments) will be varied to determine system characteristics. 

3 . 3 Computer Codes 

Analysis of the pebble bed HTGR breed/burn fuel cycle in 
the homogeneous system model described in the previous section 
required precise definition of nuclear data parameters for 
input into a code that would treat the one-dimensional diffu- 
sion-depletion problem accurately. In this section the com- 
puter codes used for this application and their respective 
input parameters are described . 

3.3.1 PDQ-7 

The PDQ-7 code [Cl] was used in the present work for the 
calculation of the operating characteristics of the breed/burn 
system model. PDQ-7 can solve the diffusion and depletion 
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equations in up to three dimensions in rectangular, cylindri- 
cal, spherical or hexagonal geometry for up to five energy 
groups. Depletion calculations and management of the neutron 
cross-section data are handled by the HARMONY system [B2], 
which provides a flexible representation of the time-dependent 
depletion equations for any nuclide chain. 

The PDQ-7 code solves the diffusion equation by discre- 
tizing the energy variable and finite differencing (central) 
the equation in space. For this application, only one-dimen- 
sional problems were treated, which the code solves by Gaussian 
elimination . 

For all calculations, PDQ-7 requires that input parameters 
define the system geometry and energy groups with respective 
chi (fission neutron) spectra, as well as the complete repre- 
sentation of nuclide chains, decay constants, initial nuclide 
number densities and microscopic neutron cross-section data. 

The problems described in Section [3.4] were all solved 
using four energy groups in an infinite slab system, with de- 
pletion parameters calculated for a static core arrangement, 
from the beginning-of-lif e to an end-of -cycle fuel burnup of 
120,000 MWD/MTHM. 

The depletion chain treatment followed 46 nuclides, for 
which the respective nuclide chains are detailed in Appendix 
[D.5], Initial nuclide number densitites were calculated as 
a function of the carbon-to-heavy metal atom ratio and the 
beginning-of-life fissile enrichment, as described in Appendix 
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[ D - 3 ] - Beginning of-life heavy metal number densities are 
given in Table [D.4] as a function of fissile enrichment, and 
all other initial concentrations are listed in Table [D.5]. 

Exposure times (in hours) used for the depletion incre- 
ments were calculated for burnup intervals as detailed in 
Appendix [D.4]. All problems were treated for burnups to 
120,000 MWD/MTHM. 

Cross-section input for PDQ-7 was provided by calculations 
using the AMPX Code System [G2], as described in Section [3.3.2], 

3.3.2 AMPX 

The AMPX Code [G2] is a modular system (NITAWL and 
XSDRNPM codes) which evaluates 123 group cross-section data 
with resonance parameters to generate coupled multi-group 
neutron cross-sections. The reference ball design used for 
number density calculations is described in Appendix [D.3]. 

The reference cell-averaged ambient temperature was 900°C [W3]. 

AMPX calculations generated all cross-section data used 
in the present work, treating the 46 nuclides listed in Table 
[D . 1 ] . For the depletion chains detailed in Table [D.7], appro- 
priate nuclide decay constants, fission product yields, energy 
releases per fission, groupwise x spectra and the resulting 
relative groupwise fluxes are given in Tables [D.8], [D.9], 

[D.10], [D.ll] and [D.12], respectively, as functions of the 
carbon- to-heavy metal atom ratio. Energy group parameters are 
given in Table [D.13], 

The four-group cross-section data generated by AMPX and 
used in PDQ-7 calculations is iisted in Appendix C. 
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3 . 4 Computational Methods for Evaluation of Homogeneous 

Zone Parameters 

As discussed in Chapter 2, evaluation of breed/burn 
operation in a homogeneous one-zone system requires computa- 
tion of system criticality and fissile enrichment as a function 
of burnup. The performance characteristics of an arbitrary 
region can be bracketed through the evaluation of several 
reference examples. 

Using the system model described in Section [3.2], the 
carbon- to-heavy metal atom ratio and the beginning-of-cycle 
fissile enrichments for the system were varied in a series of 
eigenvalue/depletion calculations. From the calculated 
information, the parameters necessary for an evaluation of 
the breed/burn system were obtained. 

For a given burn zone composition, varying the initial 
fissile enrichment will yield a series of similar curves when 
is plotted versus the fuel burnup, as illustrated in Fig. 
[3.1a]. Achieving a value of k ro > 1.0 will satisfy the first 
operating criterion: achieving system criticality. Then, a 
linear burnup model of the form: 



B E0L 

c 



a (e 



BOL 



b) 



2n 

n + 1 



(3.1) 



can be applied to the zone, where a and b are constants and n 
is the number of batch fuel loadings. For a continuously- 
fueled system, whose performance can be approached in a pebble 
bed reactor, n -»■ 00 so that: 
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FIGURE 3.1a: REPRESENTATIVE VARIATION OF 

k WITH BURNUP (linear) 

00 




FIGURE 3.1b: REPRESENTATIVE VARIATION OF 



ENRICHMENT WITH BURNUP 
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B 



E0L C = 2aU BOL - b > • 
c 



(3.2) 



The constants a and b can be calculated if batch depletion 
calculations are performed for two different initial enrich- 
ments; therefore, to achieve any design-limited end-of-life 
burnup for continuous burn zone refueling, Eq . (3.2) can be 



used to find the required z 



BOL 



Hence, for a burn zone load- 



ed with an initial fissile enrichment of e , bred in a 

JL 

breed zone, the maximum obtainable end-of-cycle burnup can be 
calculated using Eq. (3.1) for either a batch load or for 
continuous refueling . 

In its fullest development, the breed/burn cycle uses 
the heavy metal discharged from the breed zone at an enrich- 
ment £ AS y as the loaded into the burn zone at beginning- 

of-cycle, again at an enrichment e . Therefore, the value 

JL 

of s A gy charged to the burn zone must be high enough to allow 
the next cycle to achieve an equivalent breed zone discharge 
enrichment. Several successive passes through the breed zone 
may be required for a given fuel ball to achieve a desired 
enrichment, but as a very rough approximation, assuming one 
fissile nucleus is bred (and preserved) in the breed zone for 
every fissile nucleus burned in the burn zone: 



B EOL, 

bz 

B eol, 

bz 



^ - £ BOL ^ ' 
c 

^ - £ asy ^ * 



or 



(3.3) 

(3.4) 
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In evaluating this criterion, breed zone enrichments at 
breed zone end-of-cycle will be taken to be equal to the 
asymptotic values given in Table [2.5] as a function of the 
ambient spectrum. 

Burn zone enrichments can be treated more explicitly. 
Because the PDQ-7 analysis takes into account the effects of 
U235 and Pa233 production, fission product buildup, Th232 
depletion and parasitic absorption, all of which were neglect- 
ed in the simplistic analysis of Chapter 2, the actual system 
operating characteristics may differ slightly from those 
predicted previously. Concentrations of all nuclides appear- 
ing in the depletion chains are calculated by PDQ-7 at each 
time increment. Using Eq. (2.18), the enrichment of a given 
burn zone as a function of burnup can be calculated. Operating 
at a constant power density to high fuel burnups will result 
in fissile enrichments being reduced to the spectrum-averaged 
equilibrium limits at end-of-cycle, as illustrated in Fig. 
[3.1b]. If the ambient spectrum in the breed zone is much 
harder than that in the burn zone, the equilibrium enrichment 
in the burn zone will be significantly less than the asymp- 
totic limit achieved in the breed zone. Burning fuel elements 
in the burn zone to high burnups and hence to low fissile 
enrichments at end-of-life (for high carbon-to-heavy metal 
atom ratio compositions) is in keeping with the nonprolifera- 
tion objectives of the breed/burn system. Very low end-of-life 
fissile enrichments thus allow spent fuel elements to be 
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discharged directly to waste disposal, in accordance with the 
system schematic in Fig. [1.1], 

Zone compositions were evaluated using PDQ-7 in batch 
load depletion calculations, so the maximum burn zone burnups 
(achieved using continuous refueling) are twice that obtained 
by applying Eq. (3.1) (setting n = 1) to batch zone calcula- 
tions . 

In summary, the second criterion for successful operation, 
that adequate fissile material be bred in the breed zone to 
fuel the burn zone, is accomplished only if the breed zone 
end-of-life enrichment (here estimated as £ ) is: 

.0.0 JL 



'EOL. 



(%) - £ 



bz 



ASY 



(%) > B 



bz 



- EOL 



(%) 



(3.5) 



This requirement, combined with the requirement for system 

criticality are necessary (but not sufficient) conditions 

for successful breed/burn operation. The third assessment 

criterion, that of zone leakage, must be examined. 

Equation (2.39) defines the requirement that k ro >_ 1.5. 

This value applies only when the burn zone breeding ratio, 

b c , is maximized, which occurs under end-of-life asymptotic 

conditions, or when b =1.0. Since this condition (Eq. 

c 

(2.39)) defines a minimum value for k , only those combina- 

CO 

tions of breed and burn zone compositions for which the first 
two criteria for operation have been met and for which 
k >1.5 will be considered as feasible possibilities. 



BOL 
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Should all three criteria be met by any given candidate system, 
then a more detailed analysis examining system geometry in a 
two-zone heterogeneous configuration will be required. 

For this analysis, PDQ-7 depletion calculations were 
carried out for beginning-of-lif e enrichments of 2.0%, 4.5% 
and 7.0% U233, for zone carbon- to-heavy metal atom ratios of 
110, 250, 325 and 450. The analyses of k ra and e AS y as a func- 
tion of the zone burnup, previously described, were performed 
to evaluate the operating characteristics of the breed/burn 
system. Depletion calculations were run for all compositions 
through fuel burnups of 120,000 MWD/MTHM . A summary of these 
results is given in Chapter 4 . 



3 . 5 Summary 

Procedures used to analyze the characteristics of operat- 
ing breed/burn fuel cycles in pebble bed HTGRs have been pre- 
sented. A series of PDQ-7 eigenvalue/depletion runs were 
carried out on homogeneous one-zone systems to evaluate 
system performance, varying the parameters of the carbon- to- 
heavy metal atom ratio (composition) and the beginning-of-lif e 
fissile enrichment. 

Successful system operation will be achieved only if three 
necessary (but not sufficient) conditions are met: 

(a) achieving system criticality, 1.0, 

(b) achieving a fissile production rate such that 



'ASY 



(%) > B 



bz 



- EOL, 



(%) , 



accounting for zone leakage by providing adequate 
excess reactivity, k >1.5. 



and (c) 
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Satisfying all three of these objectives simultaneously will 
allow the two-zone layered pebble bed system to operate on a 
stand-alone breed/burn fuel cycle (within the limits of the 
conservative approximations used in this analysis) . 
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IV. BREED/BURN COMPUTATIONAL RESULTS 

4 . 1 Introduction 

The PDQ-7 runs made for analysis of breed/burn feasibility 
were detailed in Chapter 3. The results, to be described in 
this chapter, will characterize the operating parameters, 
consistent with core compositions and power densities as a 
function of burnup history necessary to assess the capability 
of the pebble bed HTGR system to achieve breed/burn operation. 
The analysis specifically addresses the three necessary (but 
not sufficient) criteria defined for successful breed/burn 
operation : 

(a) system criticality, >_ 1.0, 

(b) rate of fissile production, B EQL (%) £ e ASY ^ 

c bz 

and (c) accounting for zone leakage, k^ 1.5. 

4 . 2 Variation of k^ with Fuel Burnup 

Three depletion cases were treated for each homogeneous 
zone composition ( carbon- to-heavy metal atom ratios of 110, 

250, 325 and 450) , corresponding to fissile enrichments at 
the beginning-of -cycle of 2.0%, 4.5% and 7.0%. For each case, 
the depletion time increments corresponded to fuel burnups of 
2000, 6000, 30,000, 60,000 90,000 and 120,000 MWD/MTHM. Cur- 
rent materials constraints limit cycle burnups to an upper 
value of roughly 120,000 MWD/MTHM [Tl] . 
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For the twelve cases considered, the values calculated 
f or ^-oo the start of each depletion iteration are listed 

in Appendix B, Table [B.3]. As illustrated in Figs. (4.1- 
4.4] , it can be seen that for all compositions, as expected, 

^oo decreases with fuel burnup, and for each composition, high- 
er initial fissile enrichment results in higher values of 
. A step decrease in k^ is noted during the first time 
increment due to the rapid saturation of xenon-135. All 
cases exhibited a nearly linear relation between k and 

CO 

burnup beyond 6000 MWD/MTHM. 

At the beginning of life all compositions achieved criti- 
cality, with the exception of the case for which C/HM = 110 
and tgQL = 2.0%. However, all cases for which e B q L = 2.0% 
were subcritical at a burnup of only 2000 MWD/MTHM and re- 
mained so for the duration of the cycle. In these cases, a 
rise in k^ was observed at mid-cycle, apparently due to an 
increase in fissile content. This will be addressed further 
in the next section. 

The linear burnup model of Eq . (3.1) can be applied to 
each burn zone case, setting n = 1, since the calculations 
were done for batch loads. For each carbon- to-heavy metal 
atom ratio, the fissile enrichments necessary at the beginning- 
of-cycle necessary to achieve any given end-of-cycle burnup 
for either a batch-loaded or continuously-fueled system can 
then be determined. 



A linear correlation of the form: 



INFINITE-MEDIUM MULTIPLICATION FACTOR, 
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FIGURE 4.1: VERSUS FUEL BURNUP FOR C/HM = 110 
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FIGURE 4.2: VERSUS FUEL BURNUP FOR C/HM = 250 
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FIGURE 4.3: VERSUS FUEL BURNUP FOR C/HM = 325 
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FIGURE 4.4: VERSUS FUEL BURN FOR C/HM = 450 
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K co = c l B + c 0 (4.D 

was fit (linear least squares) to the linear region of each 
curve (i.e. for all burnups after initial Xel35 saturation) and 
the coefficients c q and were determined. The resulting c q 
and c 1 are listed in Appendix B, Table [B.5]. In this model, 

B is the average reactivity-limited fuel burnup. 

Equation (4.1) can be solved for the burnup at which 

A 

k = 1 . 0 , B , as : 

oo 11 



1.0 = Cl B + C Q , (4.2) 

for each burn zone burnup history shown in Figs. [4.1 - 4.4] . 
Since initial enrichments of 2.0% resulted in subcritical 
systems (for all except the first time step) for all carbon-to- 
heavy metal atom ratios examined, these cases were excluded 
from further analysis. For the remaining two burnup histories 
for each burn zone composition (corresponding to initial en- 
richments of 4.5% and 7.0%) the linear burnup relation was 
applied. Simultaneous solution of the resulting two equations 
yields the coefficients a and b for each carbon- to-heavy metal 
atom ratio. These coefficients are listed in Appendix B, 

Table [B. 6] . 

Thus, for any initial fissile enrichment, the end-of- 
cycle batch fuel burnups can be calculated; and setting n 00 
in Eq . (3.1) yields the burn zone burnup limit for continuous 

refueling. Table [4.1] lists the reactivity-limited burnups 
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TABLE 4.1 

REACTIVITY-LIMITED BURNUP* FOR 
CONTINUOUSLY-FUELED BURN ZONES 



WITH £ BOL £ ASyJ 
c bz 



BURN ZONE C/HM ATOM RATIO 
110 250 325 450 



110 

BREED 
ZONE 250 
C/HM 

ATOM 325 
RATIO 

450 



67400 97590 
28210 67390 
19920 61000 
13140 55780 



94480 86950 
65092 58500 
58880 52480 
53790 47560 



* Values shown are MWD/MTHM 

t Breed Zone e AS y values are given in Table [4.2] 



as e 



EOC 
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for continuously-fueled burn zones as a function of the breed 
zone asymptotic limits given in Table [4.2] . These asymptotic 
enrichments, being the average end-of-cycle enrichments, e , 
generated in the PDQ-7 calculations for each composition, 
approximate the asymptotic values predicted for each spectrum, 
as can be seen in Table [4.2], 

The burnups obtained vary substantially with the breed 
and burn carbon- to-heavy metal atom ratio combinations, as 
shown in Fig. [4.5] . The highest values of fuel burnups are 
obtained for a breed zone for C/HM = 110 and a burn zone 
composition in the range of C/HM = 250 to 450. 

23 

4 . 3 Variation of e with Fuel Burnup 

Having met the objective of achieving system criticality, 
the next objective is to examine the fissile production poten- 
tial. The analysis in Chapter 2 indicated that a maximum 
possible asymptotic fissile enrichment would be achieved, 
given adequate exposure time, dependent on the ambient spec- 
trum. The fissile enrichments (U233 only, neglecting U235) 
for each zone composition are listed as a function of fuel 
burnup in Table [B.4]. As illustrated in Figs. [4.6 - 4.9], 
it is clear that an asymptotic limit is approached at very 
high fuel burnups, independent of the initial zone fissile 
enrichment. A comparison is made in Table [4.2] of the average 
U233 enrichment at end-of-cycle inferred from these results 
with the asymptotic enrichment limits predicted using the 
simplistic analysis of Chapter 2. As previously discussed, 
predicted values of e ASY should correspond nearly, but not 
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TABLE 4.2 

COMPARISON OF ACTUAL END-OF-CYCLE U233 
ENRICHMENTS AND PREDICTED ASYMPTOTIC LIMITS 

C/HM ATOM RATIO 





110 


250 


325 


450 


* 

ASY 


0 - 03526 


0.02686 


0.02497 


0.02302 


'EOC 


0.03697 


0.02137 


0.01807 


0.01537 


'ASY 

'EOC 


1.048 


0.7956 


0.7237 


0.6677 



o 

c 



* Predicted using 
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FIGURE 4.5: REACTIVITY LIMITED BURNUP FOR 

asymptotic enrichment loadings 



U233 ENRICHMENT (ATOM PERCENT) 
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FIGURE 4.6 : U233 ENRICHMENT VERSUS FUEL BURNUP 

FOR C/HM = 110 



U233 ENRICHMENT (ATOM PERCENT) 
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FIGURE 4.7: U233 ENRICHMENT VERSUS FUEL BURNUP 

FOR C/HM =250 



U233 ENRICHMENT (ATOM PERCENT) 
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FUEL BURNUP (GWD/MTHM) 

FIGURE 4.8: U233 ENRICHMENT VERSUS FUEL BURNUP 

FOR C/HM = 325 



U233 ENRICHMENT (ATOM PERCENT) 
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FIGURE 4)9: U233 ENRICHMENT VERSUS FUEL BURNUP 

for C/HM =450 
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exactly, to e eol / due to the more rigorous evaluation of 
other concurrent effects by PDQ-7. 

A breed zone, initially having no fissile content, would 
eventually build up to a fissile concentration near the 
asymptotic limit, as indicated by the low-enriched (2.0% U233) 
depletion cases (which also exhibited an increase in with 
burnup due to the high fissile production rate) . It was noted 
in the previous section that an optimum combination of zone 
compositions to maximize obtainable burn zone burnup would 
use a breed zone where C/HM = 110 and a burn zone where C/HM 
= 250 to 450. Continuous refueling, which can be approached 
in pebble bed systems, would allow burn zone burnups in these 
cases approaching 100 GWD/MTHM. Recognizing that for contin- 
uous fueling very high burnups could be attained (Table [4.2]), 

the criteria for breed zone burnup, £ B qj < b E ql ' Can 

c DZ 

achieved. Additionally, it is clear from Appendix B, Table 
[B . 3] , that increasing the carbon- to-heavy metal atom ratio 
for a given initial fissile enrichment will increase k ro 

BOL 

This effect is illustrated in Fig. [4.10]. Because system 
criticality must be maintained, higher burn zone carbon-to- 
heavy metal atom ratios (v400 or greater) would therefore 
be preferred, even though the obtainable burn zone burnups 
are not quite the maximum possible. 

Breed/burn system design in the pebble bed HTGR can, 
with careful selection of zone compositions, meet the rather 
conservative criterion established here for potential fissile 
production and hence breed/burn cycle operation may be 



AT BEGINNING-OF-LIFE 
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f easible . In this analysis, the effects of poisons, control 
materials and reactor vessel materials were neglected by 
analyzing an infinite homogeneous system. Therefore, these 
results are, in fact, extremely conservative estimates, indi- 
cating that true system performance might be substantially 
poorer. The operation of a two-zone system, however, as 
described in Chapter 3, requires an evaluation which accounts 
for zone-to-zone leakage. 

4 . 4 Criticality with Zone Leakage 

Equation (2.39) required, as a minimum, that at the 
beginning-of-lif e: 



> 1.5 . (2.39) 

c 

This was necessary to account for the significant neutron 

leakage out of the burn zone (into the breed zone) . For the 

reference homogeneous systems evaluated, beginning-of-lif e 

values of k^ exceeded 1.5 only for initial enrichments of 7.0%. 

Since the asymptotic limit for fissile enrichment in a breed 

zone for which C/HM = 110 is 3.5%, as shown in Fig. [4.10], 

it can be seen that idealized core configurations yield no 

values of k _> 1.5. This indicates that it is not possible 

“bol 

to achieve a k OT high enough to accomodate zone leakage when 
using breed zone discharge fuel in the burn zone (depicted as 
the "feasible" range in Fig. [4.10]). Thus, no combination 
of the zone compositions evaluated in the present work are 
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capable of satisfying all three criteria for successful breed/ 
burn operation. For these compositions, then, it must be 
concluded that it is not feasible to operate a pebble bed 
HTGR system on the breed/burn fuel cycle. Moreover, since 
appears to attain its maximum value (for all enrichments) 
at a carbon- to-heavy metal C/HM value of 'v 4 50, it does not 
appear to be useful to consider the use of even softer spectrum 
burn zones. Referring back to Table [2.5], extrapolation to 
C/HM < 110 indicates that a harder spectrum breed zone, if 
physically practicable, might have an asymptotic enrichment 
greater than the 5% or so, which would permit k >1.5 when 

CO - — - 

cycled into the burn zone. 

4 . 5 Summary 

Analysis of the results of PDO-7 depletion calculations 
for homogeneous zones has indicated that with high fuel burnups, 
breed zone fissile enrichments approach steady state values, 
as predicted in the analysis of Chapter 2. These enrichments, 
when appropriate zone compositions are employed, will allow 
burn zones of infinite extent to achieve adequately high burn- 
ups while maintaining k ro 1.0. Considering only these cri- 
teria, breed/burn system operation appears feasible. However, 
requiring the system to overcome the effects of burn zone leak- 
age showed that none of the breed and burn zone combinations 
examined in the present work could achieve adequate values of 
k^ . Hence, the potential for successful breed/burn operation 
in pebble bed reactors in not confirmed. 
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It may be possible, however, by further reductions in the 
carbon- to-heavy metal ratio, in the breed zone, to increase 
the asymptotic enrichment limit of the breed zone (to > 5%, 
and preferably ^ 7%) and thereby make it possible to achieve 
a system 1.5. Additionally, as can be seen in Fig. [4.4], 

the use of slightly higher carbon- to-heavy metal ratios in 
the burn zone (say % 500) can produce slightly higher values 
of k ^ t and hence it may be possible to satisfy the leakage 
criterion. Reasonable compositions must be maintained, how- 
ever, because as the carbon-to-heavy metal atom ratio increases, 
the required core volume will increase for a given thermal 
output requirement. 

Another option would be to utilize the pebble bed HTGR 
in a near-breed/burn mode, by providing the additional fissile 
makeup necessary to achieve the required values of k M . Reac- 
tivity-limited burnups would thereby be increased, thus allow- 
ing the system to better utilize its potential for U233 pro- 
duction. 

The goals of non-proliferation, while best-achieved by 
a successful "stand-alone" breed/burn system, might still 
have acceptable operating parameters for a near-breed/burn 
system. Fissile U235 or U233 loaded into the system in a 
denatured form, would not violate non-proliferation standards 
and could provide the necessary reactivity to effect success- 
ful breed/burn operation in the pebble bed HTGR system. Should 
such a modified breed/burn fuel cycle be considered, however, 
the effects of plutonium production would be to be evaluated. 
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V. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 
5 . 1 Introduction 

The realization that a growing number of nonnuclear 
weapons states, desiring to provide for their own growing 
energy needs, have opted to develop civilian nuclear power 
programs, and the implied relationship between nuclear 
electric plants and weapons, has resulted in a broader inter- 
pretation of the implications of proliferation. Both national 
and international reassessments, e.g. the NASAP and INFCE 
efforts, are underway to prepare new guidelines for the devel- 
opment and application of the nuclear fuel cycle for commercial 
electrical power production. 

The transition to a breeder economy is not in prospect in 
the United States, under current governmental policy, and 
alternative options must utilize nuclear fuel in a manner so 
as to minimize the production of and/or commercial traffic in 
materials suitable for weapons applications. The present 
thesis is concerned with one such class of proliferation- 
resistant fuel cycles, designated the "breed/burn" concept. 

Fischer et al [FI] proposed the implementation of a breed/ 
burn system in a fast reactor core, using depleted uranium 
monocarbide feed fuel. After breeding sufficient fissile 
material in a hard spectrum "breed zone," fuel assemblies are 
moved to a more highly moderated region of the core where 
energy is generated by the burning of the now slightly 
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enriched fuel. 

A PPli ca tion of this concept to thermal/epithermal reactor 
systems requires the use of a composi tionally versatile system, 
capable of achieving the high fuel burnups necessary for effec- 
tive fissile conversion. 

The pebble bed HTGR offers inherent versatility in fueling 
strategies, core burnup, fuel design and neutron spectrum 
tailoring, as well as the ability to approach continuous fuel- 
ing operation by using flow-through schemes. Hence, the pre- 
sent work focuses primarily on the utilization of the U233/Th 
fuel cycle in a breed/burn configuration of the pebble bed 
HTGR. 

5 . 2 Summary of Analysis 

Fischer et al [FI] have proposed that in implementing a 
fast reactor breed/burn fuel cycle, the primary objective is 
to maximize the breeding ratio through spectrum hardening. 

More highly-moderated regions are used for fuel burning, but 
a hardened ambient spectrum in the breed zone is required, 
presumably to increase the overall system breeding ratio. 
Primarily through decreases in parasitic absorption (in 
fissile as well as structural material and fissile products) 
and increases in fertile fission rates, the breeding ratio 
does increase with hardening, but it was found in the present 
work that this is accompanied by a decrease in the asymptotic 
fissile enrichment obtained for ultra-high burnups, when the 
rate of fertile capture equals the fissile absorption rate. 
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This spectrum-dependent apparent limit is expressed by the 
relation : 



'ASY 



fertile^ fissile^ 

C cl 

1 + (a fissile^ 

C 3. 



(5.1) 



For fast reactor applications, the preference is for 
softer breed zone spectra: softer than a typical LMFBR but 
still harder than in common thermal reactors. 

Developing the breed/burn fuel cycle for thermal reactor 
configurations necessitates the use of a system capable of 
achieving regionwise spectral shifts (in order to segregate 
breed and burn zones) while being able to achieve high conver- 
sion ratios. The HTGR is the most versatile among thermal/ 
epithermal reactors currently of commercial interest. This 
inherently superior system, when employed with the U233/Th 
fuel cycle (chosen because of its outstanding neutronic 
properties in the thermal and epithermal energy ranges) 
affords the required flexibility for implementation of a breed/ 
burn cycle. Very high fuel burnups , relatively low enrichment 
requirements and the stated requirements for zoning flexibility 
made the pebble bed HTGR a logical system choice. The present 
work focused primarily on the use of this system for obtaining 
a feasible thermal/epithermal breed/burn design. 

Contrary to the situation with fast reactor breed/burn 
systems, increasing the breed region asymptotic fissile en- 
richment in the U233/Th pebble bed HTGR is accomplished by 
spectrum hardening. In addition, the values for a given 
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enrichment increase as the moderator- to-fuel ratio (carbon- 
to— heavy metal) is increased. Therefore, a system utilizing 
two zones; a very hard breed zone and a very thermal burn 
region appears most promising. 

Successful operation of any breed/burn system requires 
that: (1) system criticality be maintained, or k > 1.0, 



and (2) breed zone fissile production be adequate to fuel 

the burn zone, or, using a very rough approximation for the 

breed zone discharge enrichment, e_,. T (%) ~ (%) ~ 

tU-Li. AbY. 

bz bx 

B^ol (%) • Burn zone operation must then require an initial 
c 



enrichment, e 



gQk / no greater than that discharged from the 



'ASY ^ — B EOL ^ • A final operating 
bz c 



breed zone, or e 
criterion requires that system criticality consider the effects 
of zone-to-zone neutron leakage. In applying a leakage correc- 
tion to the macroscopic absorption cross-section, the criterion 
for criticality yields the necessity that: 



k — 1 > 1 + ^ , (5.2) 

CO — L 

a 

where e is the ’’leakage cross-section." By further manipula- 
L 

tion of Eq. (5.2), using conservative (optimistic) approxima- 
tions, one finds that: 



- 1 > 



- 1 + b 



(5.3) 



and since the burn zone breeding ratio, b , reaches a maximum 
of unity at the asymptotic limit, for successful operation. 
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accounting for zone to zone leakage / k ra must be greater than 

1.5. 



Three operating criteria have thus been established for 
successful breed/burn operation. These necessary (but not 
sufficient) conditions define the feasible range of operation 
for the system. 

To fully evaluate the practicality of pebble bed systems, 
state-of-the-art computer calculations were performed using 
the PDQ-7 code [Cl] . Four-group cross-section data was gen- 
erated for the pebble bed system using the AMPX modular code 
system (NITAWL and XSDRNPM codes) . Homogeneous one- zone 
regions were evaluated in a series of eigenvalue/depletion 
calculations, varying the zone composition (carbon-to-heavy 
metal atom ratios) and initial enrichments. Depletion runs 
were performed for fuel burnups to 120,000 MWD/MTHM . 

For the various compositions evaluated in this thesis, 
asymptotic zone enrichments are shown in the first column of 
Table [5.1]. Increasing asymptotic enrichments are observed 
for decreasing carbon-to-heavy metal atom ratios, justifying 
the preference for a hardened breed zone spectrum in this 
reactor type. Using the highest available asymptotic enrich- 



ment. 



, as the initial fissile enrichment 



C/HM = 110 

in a continuously-refueled burn zone, the burn zone end-of- 
life burnups are also shown in the first row of Table [5.1] ( 
as a function of the zone carbon-to-heavy metal atom ratio. 
It is clear that this reactor system is able to achieve high 
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TABLE 5.1 

PERFORMANCE PARAMETERS AS A FUNCTION OF 
BREED AND BURN ZONE COMPOSITION 



BURN ZONE C/HM ATOM RATIO 









110 


250 


325 


450 






b eol c ,% 


6.740 


9.759 


9.448 


8.695 






£ ASY bz ,% 




k+ 

OO 

BOL c 








110 


0.0353 


1.069 


1.246 


1.290 


1.323 


BREED 




(0.0370) 


( 1 . 110 ) 


(1.256) 


(1.286) 


(1.313) 


ZONE 


250 


0.0269 


0.971 


1.134 


1.171 


1.196 




(0.0214) 


(0.952) 


(1.095) 


(1.125) 


(1.152) 


C/HM 


325 


0.0250 


0 . 948 


1.107 


1.143 


1.166 


ATOM 




(0.0181) 


(0.910) 


(1.052) 


(1.082) 


(1.109) 


RATIO 


450 


0.0230 


0.924 


1.080 


1.113 


1.135 




(0.0154) 


(0.865) 


(1.005) 


(1.035) 


(1.062) 


* T MWD/MTHM 

ln 4 ' 

10 


using e BQL 

c 


= CftSY bz 


C/HM = 110 


for 





continuous refueling. 



For comparison, numbers in parenthesis were predicted by the 
one-group calculation of Chapter 2. 
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infinite medium fuel burnups while requiring quite low begin- 
ning-of-life enrichments. 

Values for are given in Table [5.1] as a function 

BOL 



of the combination of breed and burn zone compositions. Al- 
though most combinations result in criticality for the infinite 
burn zone (k^ > 1.0) , none of the combinations of zone composition 
reach k^ values of 1.5. Hence, based on this simplistic 
analysis, none of the compositions evaluated in this thesis 
are capable of overcoming the zone-to-zone leakage effects. 



5 . 3 Conclusions 

Successful stand-alone breed/burn operation could not 
be accomplished in U233/Th pebble bed HTGR systems for the 
ranges of compositions considered in this analysis: 

110 C/HM <_ 450. While the pebble bed system provides the 
flexibility needed to implement a successful fuel cycle, the 
effects of zone-to-zone leakage necessitates values of k^ 
that are too high to be achieved with breed zone asymptotic 
enrichments (<3.5% in the C/HM range studied). Therefore, if 
the HTGR is still to be considered as a viable breed/burn 
system, additional breed zone hardening will be necessary to 
increase the attainable fisssile enrichment: carbon- to-heavy 
metal atom ratios much less than 110 must be contemplated, so 
as to achieve asymptotic enrichments on the order or 7% or so. 
Increasing the carbon- to-heavy metal atom ratio in the burn 
zone could possibly be used to increase k^ , thereby helping 
to overcome the constraints imposed by the zone-to-zone 
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leakage. Such an action must consider, however, the additional 
reactor volume required for a given thermal output, and the 
fact that approaches a broad maximum just beyond C/HM= 450. 

For fast reactor applications, the breed/burn concept 
offers good potential for successful operation. However, 

"soft" fast reactor spectra would be most likely to meet all 
operating criteria, contrary to the proposal by Fischer et al 
[FI] . The much higher asymptotic enrichments available in 
these spectra could allow a breed/burn system to produce large 
fissile masses in the breed zone in a time span consistent 
with lower fuel burnups in the burn zone. 

A successful breed/burn system can satisfy the nonproli- 
feration objectives of minimizing the traffic in sensitive 
materials suitable for weapons application at all points in 
the fuel cycle. Additionally, the use of fertile-only makeup 
makes fuel enrichment unnecessary, while the high burn zone 
burnups attainable help reduce discharge fissile enrichments 
to levels which are insufficient to justify reprocessing. 

5 . 4 Recommendations for Future Work 

The results of this analysis showed asymptotic enrich- 
ments in fast reactor spectra were at their highest levels 
in the "soft" end of the available spectra. On the other 
hand, the highest asymptotic enrichments in the epithermal 
range were obtained in the hardest available epithermal spec- 
tra. An apparent maximum exists, therefore, for a reactor 
operating in the range between a "soft" LMFBR and a "hard" 
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HTGR. Future analyses should examine the possibility of de- 
signing a system to operate in the described range, thereby 
maximizing the potential for successful breed/burn operation. 

In addition, Karam [Kl] has proposed the use of suspended- 
bed reactors ( SBR) as alternatives to LMFBRs . These reactors 
would utilize pyrolytic carbon-coated dicarbide or metallic 
fuel kernels (i.e. HTGR fuel particles) suspended by upward 
flowing helium gas (coolant) in a cluster of pressurized tubes. 
High breeding ratios (1.2 to 1.38) are projected for this 
system, since the utilization of HTGR-type fuel kernels without 
the surrounding graphite matrix results in an extremely hard 
epithermal spectrum. Thus, it may be possible to adapt some 
of the features of the SBR to achieve a spectrum in the range 
required for maximizing the asymptotic enrichment limit. Com- 
bining the high breeding ratios and hard spectrum, and recog- 
nizing the ability to maintain composition zoning through the 
use of tube-contained fuel, this type system should be investi- 
gated carefully for its applicability to the breed/burn system 
concept . 

Less exotic HTGRs still offer exceptional potential for 
the implementation of alternate fueling schemes. Even apart 
from the SBR, conventional or pebble bed type HTGR composi- 
tions for which C/HM < 110 should be investigated, as it may 
be possible to increase the available asymptotic enrichments 
to levels high enough to produce the necessary burn zone excess 
reactivity to overcome zone- to- zone leakage effects. Only 
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after thorough analysis of all possible compositions should 
the HTGR be judged as either feasible or unfeasible for imple- 
mentation of a breed/burn fuel cycle. 

Finally, once a combination of breed/burn zones pass the 
preliminary requirements proposed in this thesis, analyses 
should proceed to the logical step of multi-zone PDQ-7 calcu- 
lations in which the zone-to-zone coupling is explicitly 
accounted for. This is essential if the extent to which zones 
can be spectrally isolated is to be established -- an essential 
prerequisite for successful breed/burn operation. 
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APPENDIX A 
ONE-GROUP CONSTANTS 

In this appendix are tabulated one-group cross-sections 
(Tables [A. 2] through [A. 5]) collapsed from the four-group 
cross-sections given in Appendix C. The four-group cross- 
sections, generated using the AMPX Code System [G2] , were 
flux weighted using the relative groupwise flux values 
listed in table [A.l]. Nuclides are identified in the cross- 
section tables by an identification number, "I#", listed 
alphabetically by isotope in Table [D.l]. 



-107- 



TABLE A . 1 

RELATIVE GROUPWISE NEUTRON FLUX* 



C/HM ATOM RATIO 



GROUP 


Au 


110 


250 


325 


450 


1 


4 . 400 


1.000 


1.000 


1.000 


1.000 


2 


5.750 


1.534 


1.533 


1.533 


1.534 


3 


5.750 


1.025 


1.194 


1.236 


1.282 


4 


5.916 


0.920 


2.305 


2.939 


3.874 



* Normalized to <J>^ for each carbon-to-heavy 



metal atom ratio. 
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TABLE A. 2: ONE-GROUP CROSS-SECTIONS FOR C/HM = 110 
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TABLE A. 4: ONE-GROUP CROSS-SECTIONS FOR C/HM = 325 
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TABLE A. 5 : ONE-GROUP CROSS-SECTIONS FOR C/HM = 450 
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APPENDIX B 

TABULATION OF DATA USED IN CHAPTERS 2 AND 4 

In this appendix, data is tabulated for the Sheaf fer 
et al [Si] cross-section correlation (Tables [B.l] and [B.2]) 
used in Chapter 2, as well as the data resulting from the 
computer runs specified in Chapter 3 and discussed in Chapter 
4, (Tables [B.3] through [B.6]). 



CONSTANTS FOR SHEAFFER [SI] CROSS-SECTION CORRELATION* 



113 - 



g 












VO 








0 












1 








rQ 






VO 






w 


CN 




i — l 


u 




CO 


OA 






OA 


O 


ro 


o 


fd 




in 


ro 






r- 


i — 1 


in 


VO 


U 


o 


i — I 


ON 


o 


o 


00 


r- 




r— 1 




o 


o 


o 


o 


o 


VD 


o 


ON 


O 

I 


g 












ro 








CD 












1 






in 


tn 






OA 






W 


00 




OA 






o 


in 






00 


CN 


r- 


VO 


X 




ON 


CO 






VO 


r- 


o 


ON 


O 


o 


rH 


1 — 1 


o 


o 




VO 


ro 


O 




o 


O 


o 


o 


o 


ON 


o 


ro 


O 



OA 

CO 

CN 


ro 


ro 

00 


O 


ro 


CN 

r- 


ro 


VO 


i — l 


in 

CN 


3 


in 


o 




■ — l 




i — l 


o 


vo 


r- 


PL| 


OA 


o 


o 


r- 


i — i 


00 


ro 


VO 


CN 




CN 


O 


ro 


i — i 


o 


1 — 1 


o 


vo 


O 



t t t 



00 




'xT 




CN 


i — l 


r- 


OA 




O 




ro 


vo 


00 


vo 




i — l 


'xT 




i — i 


O 




CN 


o 


o 


ro 


OA 


OA 


i — 1 


CN 


ro 


00 




P 


00 


o 


r- 


i — 1 




CN 


VO 


ro 


CN 






CN 


o 


CN 


o 


o 


O 


o 

1 


vo 


O 

1 
























co 






















G 


ro 




vo 






• — i 




ro 




CN 


U 


ro 


O 


00 


r- 


1 — 1 


OA 


in 


r- 




VO 


fd 


ON 


ro 


o 


00 


1 — 1 




ro 


CN 


OA 


o 


P 


P 


in 


o 


OA 


1 — 1 


ro 


CN 




i — 1 


ro 






• 


















G 




CN 


o 


CN 


CN 


O 


CN 


o 

1 


VO 


o 

I 


•H 












1 




1 




I 


CO 






















G 






















O 


ON 




vo 




OA 


OA 


00 


in 




r- 




ro 


'xT 


00 


in 


00 


in 


in 






vo 


-P 


ON 


O 


o 


00 


ro 


o 


00 


r- 


00 


00 


a 


n: 




o 


vo 


o 


00 


1 — 1 


OA 


CN 


CN 


<D 


h 




















CO 




CN 


o 


CN 


o 


o 


o 


O 


VO 


o 


1 
















1 




1 


CO 


Ph 




















CO 


w 
















U 




o 


E-i 












fd 




-P 




u 


w 












X 




K. 


U 


u 


g 






u 


1 — 1 


M-4 


1 — 1 


rd 


1 1 


-P 






:> 


UP 


O 




tn 


t) 


tn 


D 


tn 


* 



Oa 



-114- 



TABLE B. 2 

NEUTRON CROSS-SECTIONS FOR FAST SPECTRA* 



s 


cr° 2 


c 23 


2 3 

at 3 


a 28 


a 1 * 9 




3 


a 


f 


a 


a 


0.3 


0.6028 


3.739 


3.214 


0.4556 


2.623 


0.4 


0.4550 


3.306 


2.907 


0.3806 


2.401 


0.5 


0.3658 


3.005 


2.689 


0.3311 


2.242 


0.6 


0.3061 


2.780 


2.523 


0.2954 


2.121 


0.7 


0.2633 


2.603 


2.391 


0.2683 


2.023 


0.8 


0.2311 


2.459 


2.282 


0.2468 


1.941 


0.9 


0.2060 


2.338 


2.190 


0.2217 


1.872 


1.0 


0.1858 


2.235 


2.111 


0.2147 


1.813 



*From Sheaffer et al [SI] correlation, a = a^S*^. 
Cross-sections in barns. 



2.050 

1.964 

1.900 

1.849 

1.807 

1.771 

1.740 

1.713 
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TABLE B . 3 

AS A FUNCTION OF FUEL BURNUP (from PDQ-7) 



BURNUP 




C/HM ATOM RATIO 




(MWD/MTHM) 


110 


250 


325 


450 


0 


0.8863 


1.036 


1.067 


1.086 


2000 


0.8226 


0.9450 


0.9705 


0.9853 


6000 


0.8356 


0.9038 


0.9120 


0.9072 


30,000 


0.9114 


0.9091 


0.8837 


0.8335 


60,000 


0.9169 


0.8949 


0.8713 


0.8343 


90,000 


0 .9046 


0 .8724 


0.8477 


0.8102 


0 


1.164 


1.353 


1.402 


1.442 


2000 


1.103 


1.269 


1.313 


1.350 


6000 


1.081 


1.231 


1.270 


1.302 


30,000 


1.010 


1.086 


1.095 


1.090 


60,000 


0.9529 


0.9515 


0.9288 


0.8851 


90,000 


0.9180 


0.8827 


0.8523 


0.8062 


0 


1.307 


1.495 


1.548 


1.597 


2000 


1.253 


1.416 


1.463 


1.506 


6000 


1.230 


1.386 


1.431 


1.473 


30,000 


1.133 


1.257 


1.290 


1.315 


60,000 


1.036 


1.098 


1.103 


1.094 


90,000 


0.9651 


0.9601 


0.9352 


0.8893 
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TABLE B. 4 

U233 ENRICHMENT AS A FUNCTION OF FUEL BURNUP (from PDQ-7) 

BURNUP C/HM ATOM RATIO 



(MWD/MTHM) 


110 


250 


325 


450 


0 


2.0 


o 

• 

CM 


o 

• 

CM 


2.0 


2000 


1.937 


1.845 


1.828 


1.814 


6000 


2.114 


1.758 


1.681 


1.603 


30,000 


3.047 


2.019 


1.771 


1.521 


60,000 


3.422 


2.116 


1.846 


1.615 


90,000 


3.524 


2.095 


1.813 


1.574 


120,000 


3.513 


2.038 


1.753 


1.512 


0 


4.5 


4.5 


4.5 


4.5 


2000 


4.367 


4.310 


4.300 


4.292 


6000 


4.291 


4.061 


4.013 


3.964 


30,000 


4.044 


3.168 


2.944 


2.711 


60,000 


3.844 


2.487 


2.158 


1.845 


90,000 


3.720 


2.192 


1.865 


1.581 


120,000 


3.616 


2.072 


1.772 


1.524 


0 


7.0 


7.0 


7.0 


7.0 


2000 


6.840 


6.799 


6.792 


6.785 


6000 


6.667 


6.494 


6.458 


6.419 


30,000 


5.810 


5.115 


4.929 


4.729 


60,000 


4.956 


3.719 


3.388 


3.059 


90,000 


4.358 


2.772 


2.370 


1.990 


120,000 


3.962 


2.272 


1.897 


1.575 
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TABLE B . 5 

CURVE-FIT COEFFICIENTS FOR LINEAR 
k VERSUS BURNUP CORRELATION 

CO 



e BOL 

(%) 

4.5 

7.0 



COEFFICIENT C/HM ATOM RATIO 

110 250 325 

c 0 1.092 1.251 1.292 

c x -2.092E-6 -4.445E-6 -5.336E-6 

c 0 1.247 1.418 1.470 

c, -3.286E-6 -5 . 177E-6 -5.995E-6 



450 

1.326 

-6.359E-6 

1.520 

-7.009E-6 



k = c,B + c A 

OO 1 U 
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TABLE B . 6 

CURVE-FIT COEFFICIENTS FOR 
LINEAR BURNUP MODEL 



COEFFICIENT 



a 

b 



C/HM ATOM RATIO 
110 250 325 

1.256E+6 9.680E+5 9.420E+5 

1.014E-2 -1.344E-2 -1.318E-2 



450 

9.118E+5 

-1.071E-2 



B 



BOL 



= a (£ 



BOL 



- b) 



2n 



n + 1 
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APPENDIX C 

FOUR-GROUP CONSTANTS 

In this appendix are tabulated the four-group cross- 
section data generated by the AMPX code system for carbon- 
to-heavy metal atom ratio of 110, 250, 325, and 450. Refer 
to Appendix D and Section [3.3.2] for code descriptions and 
input parameters. Table [D.13] shows the neutron energy 
group structure and Table [D.l] gives the identity of each 
nuclide corresponding to the I# in the tables of this 
appendix. 
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APPENDIX D 

COMPUTER CODE INPUT PARAMETERS 
D. 1 Introduction 

In this appendix are detailed the parameters used in the 
PDQ-7 calculations described in Chapter 3. Specific spectrum- 
dependent parameters and nuclide depletion chain paramters 
were generated using the AMPX code. 

D . 2 Nuclide Identification 

Table [D.l] lists the 46 nuclides treated by the AMPX 
code system. Reference numbers are used for identification of 
nuclides in PDQ-7. 

D . 3 Number Densities 

This section describes the reference ball design used 
both in cross-section generation in the AMPX code and in 
modeling the pebble bed breed/burn system for PDQ-7 . 
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TABLE D.l 

NUCLIDE IDENTIFICATION 



NUCLIDE 


REFERENCE 


I# 5 


NUCLIDE 


REFERENCE 


I# 5 




NUMBER 4 






NUMBER 4 




Am2 4 3 


24395 


20 


Pml4 7 


14761 


12 


Carbon 


1206 


2 


Pml4 8 


14861 


14 


Csl3 3 


13355 


9 


Pml48m 


14800 


13 


Csl34 


13455 


10 


Pml4 9 


14461 


28 


Csl35 


13555 


11 


3 Pa233 


23391 


42 


Eul53 


15363 


17 


J Recvel 


1200 


6 


Eul54 


15463 


40 


Rhl03 


10345 


31 


, Eul55 


15563 


41 


Rhl05 


10545 


32 


FP 1 


10001 


22 


Sml4 9 


14962 


15 


FP 2 


10002 


21 


Sml50 


15062 


38 


Helium 


402 


1 


Sml51 


15162 


16 


1135 


13553 


4 


Sml52 


15262 


39 


Kr82 


8236 


29 


Silicon 


1428 


7 


Kr 83 


8336 


30 


Agio 9 


10947 


33 


Ndl43 


14360 


36 


Th2 32 


23290 


23 


Ndl45 


14560 


37 


U233 


23392 


24 


Ndl4 7 


14760 


3 


U234 


23492 


25 


Np239 


23993 


27 


U235 


23592 


26 


Oxygen 


1608 


5 


U236 


23692 


18 


Pu239 


23994 


44 


U238 


23892 


43 


Pu240 


24094 


45 


Xel31 


13154 


34 


Pu241 


24194 


46 


Xel33 


13354 


35 


Pu2 4 2 


24294 


19 


Xel35 


13554 


8 



Slov/ly-saturating fission products (lumped) 

2 

Non-saturating fission products (lumped) 

3 

Reciprocal-velocity correction to account for graphite 
impurities [W3] 

4 

Combined mass number and atomic number for identification 
in PDQ-7 

5 

Identification number for AMPX code and cross-section 
tables of Appendices [A] and [C] 



- 150 - 



D.3.1 Nomenclature 



C/HM = Ratio of carbon number density to heavy metal 
number density 

e = Fissile enrichment 
N 2 3 

N 2 3 + N° 2 



N 



A 



1 



M 



1 



P.F. 



V 

V. 



1 



V 

x 



1 



V’ 

x 



2 



N 



x 



2 



N 



x 



Avogadro's number 



0.6022 



cm - atoms 
b-mole 



Atomic mass of indicated nuclide, [gm/mole] 
Density of indicated nuclide, [gm/cm 3 ] 

Packing factor of balls 

Volume of balls = o 61 
Volume of reactor 

3 

Volume of single ball, [cm ] 

Effective volume of ball (cell-averaged) , [cm 3 ] 

Volume of indicated nuclide in ball, [cm 3 ] 

Effective volume of nuclide in ball (cell- 
averaged) , [cm 3 ] 

Atom density of indicated nuclide in ball, 

[b 1 cm *] 



Effective atom density of indicated nuclide 
in ball [b ’em -1 ] 



= Radius of ball, [cm] 



^Subscript T, C or U indicates Th 232 0 2 , graphite or U 233 0 2 , 
respectively . 

2 

Subscript T, C, U or 0 indicates Th232, graphite, U233 or 
oxygen, respectively. 
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D . 3 . 2 Material Properties 

Calculation of initial number densities, as a function of 
U233 enrichment, involved considerations of the volume fractions 
of the particular materials within a fuel cell. Fuel cells 
were defined as the volume of one fuel ball divided by the 
packing factor, P.F. Properties of the major substances in 
the fuel balls are shown in Table [D.2] 

The makeup of the fuel balls is illustrated in Figure [D.l], 
appropriately dimensioned. Packing fractions of fueled grains 
in the meat are listed in Table [D.3]. Number densities for 
all materials present at beginning-of-lif e are given, and 
all others appearing in the depletion chains are initially 



zero. 



TABLE D.2 



BALL MATERIAL 


PROPERTIES [W3] 




Material 


p 


M 




F x 


X 




(gm/cm 3 ) 


(gm/mol 


Heavy metal oxides 


9.5 


264 


Graphite : 






grain coating 


1.5 


12.0 


meat 


1.0 


12.0 



outer shell 



1.75 



12.0 
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FIGURE D.l: SECTIONAL VIEW OF FUEL BALL 

(not to scale) 
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TABLE D . 3 

PACKING FRACTIONS OF GRAINS IN MEAT [W3] 

C/HM PACKING 

ATOM RATIO FRACTION* 



110 


0.39073 


250 


0.168 


325 


0.127 


450 


0.091987 



*Packing fraction - 


(Number of grains) (Volume per grain) 
Volume of Meat 
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D-3.3 Calculation of Number Densities 
The volume of a fuel element is given by: 

V B = I 771 * 3 = 113.1 cm 3 , 

so the effective volume (cell-averaged) is: 

V T 



(D.l) 



VI = 



B =185.4 cm 3 . 



(D.2) 



B P.F, 

Number densities for thorium, oxygen, uranium, and carbon were 
determined using the following relationships: 



N 



h HM a 



HM 



M 



HM 



p c n a 

N = ----- 
C M„ 



(D. 3) 



(D. 4 



= 



P 0 N A 



0 M 



(D.5) 



0 



where the densities, p, are "smeared" throughout the ball. 
In addition, 



N. 



N, 



C/HM = 



N + N N 
T lN U HM 



(D.6) 



and 



e = 



N 



U 



= N U 



( D . 7 ) 



N + N N 
T U HM 



To find the smeared densities for heavy metal and carbon: 






= P 



HM 



HM HM0 2 \M hm0 ^ 



V 



Kernel 



\ 



V 



Grain 



V Grain | f^Meat 



, V 



Meat j \ Ball 



(D. 8) 
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p = p • II - VKern el U V Grain | 
P C p Grain 1 1 V„ Hv 1 



\ / 



+ p Meat 1 V 



Grain 

V-. . \ 

Grain 



Meat 



V 



Meat 

\ 



V 



Meat I 



V 



Ball i 



+ P Outer Shell l 1 V 



Meat 
V Ball 

V. \ 



Meat 

Ball 



I ■ 



(D.9) 



Since N ' = N (P.F.), the effective number densities (cell- 
averaged) can be summarized as follows: 



V = 



2 M. 



Iv„ ^ \ N,. 

Kernel I (p>F>) A 



HMO. 



V 



Ball i 



M. 



HM 



(D.10) 



n t’ 



= (1 - e) P 



M HM \/ V Kernel) 



N, 



HM0 2 M HMO. 



V 



(P.F.) 



Ball 



M 



HM 



(D . 11) 



N, 



N c ' 



N 



(P.F.) 



HM 





V Kernel\ 


| V Grain | 


P Grain 


i V Grain ) 


\ V Ball ) 



V 



+ p Meat I ^ 



Grain 



V 



Meat 
V 

+ p Outer Shelll 1 “ V 



\| 


fv Meat\ 




y V Ball / 



Meat 



Ball 



(D . 12 ) 



In addition to the above, initial number densities were 
used as input for silicon; the recvel correction (reciprocal 
velocity [W3]) for graphite impurity is used at a concentration 
of approximately 2 ppb carbon, so that the effective number 
densities (cell-averaged) become: 



N 



recvel 



(2.05E-9)N C * 



(D. 13) 



N ' 

helium 



8 . 308E-5 



(D . 14 ) 
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N ' . , . 
silicon 



(9. 14E-3) 



^Meat | / V Grain 

V V 

Ball l\ Meat. 



Representative values for N y ' and N^' are given in 
Table [D.4], while values for N^ M , Nq , N£ , N ^. ecvel / 

Silicon' and N helium are listed in Table [D.5], 



.15) 
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TABLE D . 4 

INITIAL CELL-AVERAGED HEAVY METAL NUMBER DENSITIES* 



C/HM 

ATOM RATIO 


£(%) 


U233 


Th232 


110 


2.0 


7.9408E-6 


3 . 8910E-4 


110 


4.5 


1 . 7867E-5 


3. 7918E-4 


110 


7.0 


2. 7793E-5 


3. 6925E-4 


250 


2.0 


3. 4143E-6 


1. 6730E-4 


250 


4.5 


7. 6821E-6 


1. 6303E-4 


250 


7.0 


1 . 1950E-5 


1. 5876E-4 


325 


2.0 


2 . 5810E-6 


1 . 2647E-4 


325 


4 . 5 


5. 8073E-6 


1 . 2324E-4 


325 


7.0 


9. 0336E-6 


1. 2002E-4 


450 


2.0 


1.8695E-6 


9.1604E-5 


450 


4.5 


4 . 2063E-6 


8.9267E-5 


450 


7.0 


6.5431E-6 


8 . 6930E-5 



* t ,_-l -1 

* In b cm 



INITIAL CELL NUMBER DENSITIES* 
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D . 4 Exposure Times 

This section describes the calculation of exposure times 
(hours) used in PDQ-7 depletion calculations for fuel burnups 
up to 120,000 MWD/MTHM. Calculated exposure times are tabu- 
lated as a function of the carbon- to-heavy metal atom ratio 
in Table [D. 6] . 

D.4.1 Nomenclature 

q' ' 1 = Average core power density 

3 

= 5.0 watts/cm 

N HM = Cell-averaged heavy metal atom density, 
[b -1 cm -1 ] 

t = Exposure time at full power, [hours] 

N^ = Avogadro ' s number 

2 

= 0.6022 cm -atoms/b-mole 

M r „, = Atomic mass of heavy metal nuclides 
HM 

- 232 gm/mole 

B = Fuel burnup, [MWD/MTHM] 

3 

V„ = Volume of relevant reactor portion, [cm ] 
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D . 4 . 2 Calculation of Exposure Times 

Exposure times, t [hours], were calculated to correspond 
to specific values of fuel burnup, B [MWD/MTHM] . As a function 
of the cell-averaged heavy metal number density, [b ‘'"cm ■*"] , 

the fuel burnup is given by: 



q' ' *N A t 

B = (2.4E-23) MWD/MTHM 

n hm m hm 



(D . 16) 



so that the exposure time, t [hours], can be expressed as 

t [hours] = (1.8492E+3)B [MWD/MTHM] N^ M [b _1 cm _1 ] (D.17) 

where can be found by summing Equations (D.10) and (D.ll). 

Representative values for the exposure times are given in 
Table [D.6] as a function of C/HM atom ratio. 
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TABLE D . 6 
EXPOSURE TIMES* 

C/HM ATOM RATIO 



B (MWD/MTHM) 


110 


250 


325 


450 


100 


7. 3421E+1 


3. 1568E+1 


2.3864E+1 


1 . 7 285E+1 


1000 


7. 3421E+2 


3 . 1568E+2 


2 . 3864E+2 


1. 7285E+2 


2000 


1. 4684E+3 


6. 3136E+2 


4 . 7728E+2 


3.457 OE+2 


3000 


2. 2026E+3 


9.4704E+2 


7 . 1592E+2 


5 . 1855E+2 


6000 


4. 4053E+e 


1. 8941E+4 


1. 4318E+3 


1. 0371E+3 


10,000 


7 . 3421E+3 


3 . 1568E+3 


2 . 3864E+3 


1. 7285E+3 


15,000 


1.1013E+4 


4. 7352E+3 


3 . 5796E+3 


2 . 5927E+3 


20,000 


1. 4684E+4 


6. 3136E+3 


4 . 7728E+3 


3 . 457 0E+3 


25,000 


1. 8355E+4 


7 . 8920E+3 


5 . 9660E+3 


4 . 3212E+3 


30,000 


2. 2026E+4 


9 . 4704E+3 


7 . 1592E+3 


5 . 1855E+3 


45,000 


3. 3040E+4 


1. 4206E+4 


1. 0739E+4 


7.7783E+3 


60,000 


4. 4053E+4 


1. 8941E+4 


1 . 4318E+4 


1 . 0370E+4 


75,000 


5. 5966E+4 


2. 3676E+4 


1. 7898E+4 


1.2964E+4 


90,000 


6 . 6079E+4 


2. 8411E+4 


2 . 1478E+4 


1 . 5556E+4 


100,000 


7. 3421E+4 


3.1568E+4 


2. 3864E+4 


1.7285E+4 


120,000 


8 . 8106E+4 


3. 7882E+4 


2 . 8 637E+4 


2. 0742E+4 



*Times in hours. 

For comparison, note that 8766 hours = 1 year. 
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D.5 Depletion Treatmen t 

D.5.1 Depletion Chains 

Depletion chains for the 46 nuclides considered are shown 
in Table [D. 7], PDQ-7 solves the sequence of chains in order. 
Chain 1 through Chain 19, for each pass through the depletion 
mesh. In Table [D.7] Chains 1 through 14 are depleting. 
Underscoring indicates the existence of a fission yield, given 
in Table [D.9], for the indicated nuclide. Reactions are in- 
dicated as follows: 

(1) : (B~) 

(2) : (n, y) 

(3) : (n, y) plus (8~) 

D.5.2 Nuclide Decay Constants 

In Table [D.8] are tabulated all nuclide decay constants 
used in the subject PDQ-7 calculations. 

D . 5 . 3 Fission Product Yields 

Table [D.9] lists the yields of each fission product 
treated for six fertile and fissile nuclides. 

D.5.4 Energy Release per Fission 

The energy release per fission for all heavy metals treated 
by the AMPX code are given in Table [D.10] . 

D . 5 . 5 Chi (Neutron Yield) Spectra 

As a function of the carbon-to-heavy metal atom ratio. 



the groupwise distributions for fission neutrons, X/ are given 
in Table [D. 11] . 
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TABLE D.7 (page 1 of 4 ) 
DEPLETION CHAINS 



CHAIN 1: U23 3 Chain 




Pa233 

U) \ ( 2 ) ( 2 ) ( 2 ) 

U233 >• U234 > U235 >U236 



CHAIN 2: Xel35 Chain 





1135 


( 1 ) 


1 




Xel35 


( 1 ) 


4 



Csl35 

CHAIN 3: Sml49 Chain 1 



Ndl47 

( 1 ) \ ( 2 ) ( 2 ) 

Pml47 Pml4 8m — =>■ Pml49 

(iT 



( 2 ) 



->■ Sml52 



( 3 ) 



( 2 ) ( 
Sml49 > Sml50 - 



( 2 ) ( 2 ) 

Eul53 Eul54 >■ Eul55 



2 ) 

— >• Sml51 
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TABLE D. 7 (page 2 of 4 ) 
CHAIN 4: Sml49 Chain 2 



( 2 ) 

— > Sml51 



CHAIN 5 : Kr83 Chain 

( 2 ) 

Kr82 >Kr83 

CHAIN 6 : Rhl03 Chain 

Rhl03 

( 1 ) 

V 

CHAIN 7 ; Rhl05 Chain 

Rhl0 5 

( 1 ) 

V 

CHAIN 8 : Agl09 Chain 

Agl09 

( 1 ) 

V 



Ndl47 



( 1 ) 



( 2 ) ( 2 ) 

Pml47 > Pml4 8 > Pml49 

or- 



( 2 ) 

■Sml49 > Sml50 



( 2 ) 

— >Sml52 



(3) 



( 2 ) ( 2 ) 

k Eul53 >Eul54 ^Eul55 
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TABLE D. 7 (page 3 of 4 ) 



CHAIN 9: Xel31 Chain 



Xel31 



(i: 



CHAIN 10 : Fission Product Chain (lumped) 

( 2 ) 

FP1 5- FP2 



CHAIN 11: Ndl43 Chain 



Ndl4 3 



( 1 ) 



CHAIN 12 : Xel33 Chain 

Xel33 

(1) 1 ( 2 ) 

Csl33 >* Csl34 

CHAIN 13: Ndl45 Chain 



Ndl45 



( 1 ) 



CHAIN 14: Pu239 Chain 



U238 

( 3 )^ 



Np239 



( 1 ) 



I (2) (2) (2) 

Pu239 ^ Pu240 ^ Pu241 ^Pu24 2 



(3 






Am243 
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TABLE D.7 (page 4 of 4) 


CHAIN 15: 


Helium Chain (non-depleting) 



Helium 



CHAIN 16: 


Carbon Chain (non-depleting) 



Carbon 



CHAIN 17: 


Oxygen Chain (non-depleting) 



Oxygen 



CHAIN 18: 


Recvel Chain (non-depleting) 



Recvel 



CHAIN 19: 


Silicon Chain (non-depleting) 




Silicon 
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TABLE D . 8 (page 1 of 2) 
NUCLIDE DECAY CONSTANTS 



CHAIN 

# 


I# 


REFERENCE 

NUMBER 


NUCLIDE 


1 


23 


23290 


Th232 


1 


42 


23391 


Pa233 


1 


24 


23392 


U233 


1 . 


25 


23492 


U234 


1 


25 


23592 


U235 


1 


18 


23692 


U236 


2 


4 


13553 


1135 


2 


8 


13554 


Xel35 


2 


11 


13555 


Csl35 


3,4 


3 


14760 


Ndl47 


3,4 


12 


14761 


Pml47 


3 


13 


14800 


Pml48m 


4 


14 


14861 


Pml4 8 


3,4 


28 


14961 


Pml4 9 


3,4 


15 


14962 


Sml49 


3,4 


38 


15062 


Sml50 


3,4 


16 


15162 


Sml51 


3,4 


39 


15262 


Sml52 


3,4 


17 


15363 


Eul53 


3,4 


40 


15463 


Eul54 


3,4 


41 


15563 


Eul55 


5 


29 


8236 


Kr8 2 


5 


30 


8336 


Kr83 


6 


31 


10345 


Rhl03 


7 


32 


10545 


Rhl05 


8 


33 


10947 


Agl09 



*X, [sec 1 I 

1. 58E-18 

2 . 93E-7 

1 . 358E-13 

8. 856E-14 

3.08E-17 

9 . 189E-16 

2.87E-5 

2.09E-5 

1 . 098E-14 

7.32E-7 

8.29E-9 

1.98E-7 

1.49E-6 

3 . 63E-6 

5. 491E-23 

Stable 

2.75E-10 

Stable 

Stable 

1. 37E-9 

1. 16E- 8 

Stable 

Stable 

Stable 

5.36E-5 

Stable 



TABLE D . 8 (page 2 of 2) 



9 


34 


13154 


10 


22 


10001 


10 


21 


10002 


11 


36 


14360 


12 


35 


13354 


12 


9 


13355 


12 


10 


13455 


I 3 


37 


14560 


14 


43 


23892 


14 


27 


23993 


14 


44 


23994 


14 


45 


24094 


14 


46 


24194 


14 


19 


24294 


14 


20 


24395 


15 


1 


402 


16 


2 


1206 


17 


5 


1608 


18 


6 


1200 


19 


7 


1428 



For stable nuclides, X 



Stable 
Stable 
Stable 
Stable 
1. 52E-6 
Stable 
9.55E-9 
Stable 
4. 87E-18 
3. 41E-6 
9 . 016E-13 
3. 249E-12 
1.68E-9 
5. 795E-14 
2. 871E-12 
Stable 
Stable 
Stable 
Stable 
Stable 

0.1E-20 sec 



Xel31 

FPl 

FP2 

Ndl43 

Xel33 

Csl3 3 

Csl34 

Ndl4 5 

U238 

Np239 

Pu249 

Pu240 

Pu241 

Pu24 2 

Am243 

Helium 

Carbon 

Oxygen 

Recvel 

Silicon 



was used. 



I# 

4 

8 

28 

36 

22 

21 

29 

30 

31 

32 

33 

34 

35 

37 

17 

41 

3 

16 

39 
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TABLE D . 9 

FISSION PRODUCT YIELDS 



Th232 


U233 


FISSION 

U235 


0.0539 


0.0562 


0.0617 


0.0021 


0.0022 


0.0024 


0.00945 


0.0077 


0.0113 


0.08 


0.0591 


0.0603 


0 . 301 


0.309 


0.298 


1.137 


1.119 


1.26 


0.01 


0.0069 


0.0028 


0.0199 


0.0114 


0.00544 


0.0016 


0.016 


0.03 


0.0005 


0.005 


0.009 


0.00055 


0.0004 


0.0003 


0.017 


0.0341 


0.0293 


0.037 


0.0588 


0.659 


0.066 


0.0338 


0.0398 


0.00038 


0.0015 


0.00169 


0.00009 


0.0003 


0.00033 


0.038 


0.0193 


0.0236 


0.002 


0.0035 


0.0044 


0 . 00085 


0.0022" 


0.00281 



SOURCE 






U238 


Pu2 39 


Pu241 


0.0578 


0.0693 


0.0626 


0.0022 


0.0027 


0.0024 


0.021 


0.013 


0.012 


0.48 


0.046 


0.056 


0.273 


0.324 


0.324 


1.426 


1.456 


1.456 


0.0005 


0.0017 


0.0007 


0.0009 


0.0029 


0.001 


0.062 


0.056 


0.062 


0.037 


0.055 


0.058 


0.003 


0.0155 


0.031 


0.023 


0.038 


0.034 


0.052 


0.069 


0.06 


0.039 


0.031 


0.036 


0.0038 


0.0034 


0.002 


0.0013 


0.0016 


0.001 


0.028 


0.0205 


0.022 


0.009 


0.008 


0 . 005 


0.006 


0.006 


0.003 
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TABLE D.10 

ENERGY RELEASE PER FISSION 

FISSIONABLE K 

NUCLIDE [watt-sec/fission] 



Th2 3 2 


3.18 


X 


10 -1 1 


Pa233 


3.19 


X 


10" 1 ’ 


U23 3 


3.19 


X 


10 _1 1 


U234 


3.19 


X 


10" 1 1 


U2 3 5 


3.21 


X 


lCT 1 1 


U236 


3.26 


X 


10 _1 1 


U238 


3.31 


X 


lCT 1 1 


Np239 


3.27 


X 


10 _1 1 


Pu239 


3.34 


X 


io " 1 1 


Pu24 0 


3.36 


X 


lCT 11 


Pu24 1 


3.37 


X 


IO - ' 1 


Pu24 2 


3.37 


X 


10"’ ! 


Am24 3 


3.39 


X 


10" 1 1 



TABLE D.ll 



CHI 


(NEUTRON YIELD) 


SPECTRA 






C/HM ATOM 


RATIO 




110 


250 


325 


450 


9 . 66727E-1 


9 . 66727E-1 


9. 66727E-1 


9. 66727E-1 


3. 32743E-2 


3.32745E-2 


3. 32746E-2 


3. 32743E-2 


4 . 90742E-9 


5. 67340E-9 


5 . 91979E-9 


5 . 17704E-9 



4 . 98079E-14 



5 . 75821E-14 



6 . 00830E-14 



5 . 25449E-14 
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D.5.6 Neutron Flux 

Table [D.12] tabulates the groupwise relative neutron 
flux generated for each carbon-to-heavy metal atom ratio by 
the AMPX code. Group parameters are given in Table [D.13]. 
These flux values were used to collapse the four-group con- 
stants given in Appendix C to one-group, for use in the 
analysis of Chapter 2. 
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TABLE D . 12 

GROUPWISE RELATIVE NEUTRON FLUX* 



C/HM ATOM RATIO 



FLUX 


110 


250 


325 


450 


*1 


1.000 


1.000 


1.000 


1.000 


(j> 2 


1.534 


1.533 


1.533 


1.534 


(p 2 


1.025 


1.194 


1.236 


1.282 




0.920 


2.305 


2.939 


3.874 



* Normalized by using <f>g/<f>g for each C/HM 
atom ratio. 
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Upper Energy (eV) 
Lower Energy (eV) 
Lethargy width, Au 



TABLE D.13 
ENERGY GROUPS 



Group 1 

1.4918E+7 
1. 8316E+5 
4.400 



Group 2 

1 . 8316E+5 
5 .8295E+2 
5.750 



Group 3 

5.8295E+2 

1.8554E+0 

5.750 



Group 4 

1 . 8 554E+0 

5.00E-3 

5.916 
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